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DIGITAL COMPUTER ANALYSIS OF CONDENSATION 
IN HIGHLY EXPANDED FLOWS 



SUMMARY 

The application of the IBM 7090 Digital Computer to the theoretical 
prediction of condensation in highly expanded flows is presented. The 
equations of the spontaneous nucleation theory of Frenkel (13) are com- 
bined with the steady one -dimensional diabatic flow equations for a solu- 
tion of the expansion of a pure vapor. A digital computer program for 
the solution of these equations is compiled and presented. 

The theoretical prediction of the condensation of nitrogen is com- 
pared to experimental results and variations in specific heat, latent 
heat of vaporization, surface tension, and rate of expansion are inves- 
tigated. The theoretical calculation is applied to metal vapors and 
the results for copper and zinc vapors are presented. 

The results indicate that the theoretical solution gives a reason 1 - - 
able prediction of the condensation in highly expanded flows. The 
degree of supersaturation increases with an increased rate of expan- 
sion and for a proper set of initial conditions "condensation free" flow 
is obtainable. The rate of expansion and the surface tension are the 
most critical parameters in the equations for condensing flow. Vari- 
ations of specific heat and latent heat of vaporization show only minor 
effects on the end result. 
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I. INTRODUCTION AND SURVEY OF THE LITERATURE 



The phenomenon of the spontaneous condensation of a gas or vapor 
into the liquid phase has been the subject, of numerous theoretical and 
experimental investigations for the past three-quarters of a century. 
The condensation of an expanding vapor was first noted in 1887 by 
R. von Helmholtz \ however, this phenomenon did not enter the field of 
aeronautics until about 1940 when its occurrence was noted in the 
supersonic wind tunnel. 

In a supersonic nozzle, the gas or vapor undergoes an isentropic 
expansion which results accordingly in a decrease in the temperature 
and pressure. Generally the saturation vapor pressure decreases with 
temperature more rapidly than does the static pressure of the isen- 
tropic expansion, hence saturation conditions are approached. If the 
rate of expansion is very rapid the flow will pass through the point at 
which saturation temperature and pressure are reached without any 
condensate being formed. The vapor is then in a non-equilibrium state, 
supersaturated, in which the existing vapor pressure is higher than the 
corresponding equilibrium saturation pressure. 

The first person to note that a supersaturated vapor can exist was 
R, von Helmholtz, who showed that saturated steam expanding through 
an orifice into the atmosphere would remain clear for some distance 
before becoming cloudy, Stodola (1) furthered this by devising experi- 
ments on superheated steam which would detect the onset of condensa- 
tion and provide evidence of super saturation, 
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In the late 1940's and the early 1950's an intensive series of inves- 
tigations were conducted to determine the condensation of air, water 
vapor in air, and nitrogen. Many of these were conducted at the Cal- 
ifornia Institute of Technology, and worthy of note are the investiga- 
tions of Head (2) on water vapor in air, Willmarth and Nagamatsu (3) 
on nitrogen, and a summary report of GALCIT work by Nagamatsu (4). 

A series of investigations of air in high speed wind tunnels conducted 
during this era include the works of Stever and Rathbun (5), Wegener 
and Reed (6), McLellan (7) and Bogdonoff and Lees (8). 

There was a lag of interest in this area in the middle to late 1950’s. 
However, with the advent of the nuclear space propulsion system and 
the use of metallic vapors as the working fluid for hypersonic wind 
tunnels, the interest has been renewed. The Office of Naval Research 
sponsored the recent review of existing theories by Courtney (9) in an 
effort to stimulate this interest. Many of these theories are built 
around the spontaneous self-nucleation concept, and they base the 
collapse of the supersaturated state upon the determination of the rate 
of droplet growth to the critical drop size. The complexity of the 
nucleation phenomenon requires that many simplifying assumptions 
be made in order to evaluate the problem by hand calculations. 

Courtney (9) attempts to obtain a more complete evaluation of the 
nucleation theories by carrying out the calculations on a digital 
computer. 

A recent study of air in hypersonic wind tunnels by F. L. Daum (10) 
indicates that, as yet, there is no proven available theory for predict- 
ing the amount of supersaturation expected in a hypersonic nozzle from 
which may be determined the point of onset of condensation. 
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Until recently none of the investigations had included the problem 
of metal vapors. The investigation by Hill, et al, (11) in 1962 adapted 
the existing theory to the case of metal vapors and provided results 
which indicate that this is a valid extension of the theory. The devel- 
opment of the nucleation equations and flow equations did not prohibit 
their application to metals, but since the need had not existed the 
application was not attempted. The authors noted have not mentioned 
any adaptation of the digital computer to the hypersonic flow problem. 
Hill suggests that this is the next logical step and that a computer pro- 
gram should be used to investigate the effect of variations in surface 
tension, nozzle length, area distribution, inlet conditions, and other 
parameters as well. 

Purpose of the Present Investigation 

At the inception of this study there was no known application of 
condensation theory to the case of metal vapors. The Hypersonic 
Wind Tunnel Group of The University of Michigan was undertaking 
the experimental investigation of the condensation of metal vapors in 
highly expanded flows, and thus arose the need of a method for the 
theoretical prediction of the onset and magnitude of condensing flow. 

Of greater importance was the need for an understanding of the 
effect a variation in any of the various parameters might have on the 
end result. Thereby arose the need for a digital computer program 
from which numerous results can be obtained for multiple combina- 
tions of the various parameters. 

The extent of the present study is to: 

(a) review the literature to determine which of the more promis- 
ing theories will most readily adapt to the IBM 7090 digital 
computer, 
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(b) write the program and test against some known experimental 
results to validate the method, 

(c) obtain theoretical results for several metal vapors under 
reasonable hypersonic test conditions, 

(d) vary surface tension, specific heat, heat of vaporization, and 
nozzle parameters to show the magnitude of influence on the 
problem. 

The first approach was one toward a general problem with all 
quantities variable and with multicomponent working fluids. This 
rapidly exceeded the capabilities of the computer, and of the author, 
and dictated the simplified approach which is employed. The system 
will be developed for one-dimensional expansion of pure vapors which 
obey the perfect gas law and which have a constant ratio of specific 
heats. Initially the specific heat, latent heat, and surface tension will 
be taken as constant and later will be varied to observe the effect. 



II. THERMODYNAMIC FUNDAMENTALS 



A. Isentropic Expansion of a Vap or 

The assumption that the vapor follows the perfect gas law is, in 
general, well justified and the pressure and temperature can be related 
by 



P 



P 



o 




( 2 . 1 ) 



Any vapor whose latent heat is large, i. e. , L > C T, will approach 

P 

saturation as the temperature is decreased by an isentropic expansion, 
and if the expansion is rapid enough there will be a thermal lag so that 
a supersaturated region will be reached. 

The Clausius-Clapeyron equation can be used to relate the variation 
of the saturation vapor pressure with temperature to the latent heat of 
vaporization. 
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(2. 2) 



where p^ is the saturation vapor pressure corresponding to a partic- 
ular temperature and for an infinite droplet radius, L(T) is the latent 
heat of vaporization, and v^ and v^ are the specific volumes of the 
vapor phase and liquid phase, respectively. 

Since the liquid occupies a very small portion of the total volume, 

we will assume v - v T = v and v = — . Assuming L = constant, 
v L v v p 

1 V 

equation (2. 2) becomes 
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where, p. - the molecular weight of the vapor 
R -• universal gas constant. 

In order to compare the slope of the isentrope with (2. 3) we will 
differentiate (2. 1) logarithmically f 
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from which, 
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dT 
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(2. 5) 



by virtue of v - C /C and R/ u - C - C . The ratio of slopes is 
J p v p v 

therefore 



dp , 

_ oo i dp 

dT ! dT 



C T 
P 



( 2 . 6 ) 



Thus for the case of L/T > the saturation vapor pressure decreases 
with decreasing temperature more rapidly than does the vapor pressure 



in an isentropic expansion, and the supersaturation ratio. 



00 ' 



contin- 



ually increases until saturation and, eventually, super saturation are 
obtained. 

Head (2) gives slope ratios for several vapors and shows a typical 
value of 2 for nitrogen nozzles operating at normal reservoir conditions. 
Similarly, he shows that if a vapor had a small latent heat of vaporiza- 
tion, L/T < Cp, an isentropic compression is necessary for the approach 
to saturation conditions. 
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Upon reaching the saturation point the fluid has a choice of two 
routes by which to continue the expansion. One, the saturated equilib- 
rium expansion, involves a gradual process of condensation which com- 
mences at the saturation point and maintains thermodynamic equilibrium 
throughout the expansion. The other, a supersaturated expansion, is 
characterized by a delay in the onset of condensation to some point down- 
stream of the saturation point at which the vapor temperature is well 
below the corresponding saturation temperature. 

Before continuing further with these alternatives it will be neces- 
sary to develop the diabatic flow equations which account for the effect 
of the condensate. 

B. One -Dimensional Steady Flow Equations With Condensation 

The equations have been derived by many authors; however, we are 
interested here in a specific case which, includes the following assump- 
tions: 

1. The vapor phase is a pure vapor and may be treated as a per- 
fect gas. 

2. The condensed mass is a liquid, is uniformly distributed 
throughout the gaseous components, and has the same speed 
and temperature as the stream. 

3. The volume of the condensed phase is negligible compared to 
the total volume. 

4. The nozzle flow is frictionless and is without heat transfer. 

5. The saturation curve for the vapor may be adequately approxi- 
mated by the Clausius - Clapeyron equation. 

6. The latent heat of vaporization and the specific heat at constant 
pressure remain constant. 
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Assumption 2 requires some discussion, in that the crystallization 
of the liquid drops or the sublimation directly to the solid state is prob- 
able at the lower temperatures. In the area concerned, the available 
experimental data on the properties of the vapor and the liquid phases 
is extremely limited, and for the crystalline phase it is practically 
non-existent. Highly extrapolated curves are the only sources, in many 
instances, so the liquid assumption will be as good as the available 
data. The zero drag approximation is supported by Stever (12) in cit- 
ing rumerous investigations which verify high fractional condensation 
rates. From the equations of Frenkel (13) it is noted that a small crit- 
ical drop radius is necessary for a high condensation rate. In high 
speed flow the drops are in the nozzle for only a short period of time, 
so there is only time for a small increase in size due to growth. Thus, 
for high speed expansions with supersaturation, the drop size will be 
small, and the logical assumption is that the condensate and the vapor 
travel at the same velocity. 

Only the major steps in the development of the flow equations for 
a pure vapor are presented in this section. A more complete develop- 
ment is presented in Wegener (14). 

The continuity equation can be written in its usual form, since the 
volume of the liquid is negligible and the vapor and drops travel at the 
same velocity. With p taken as the density of the mixture and m as the 
total mass which passes a given location in unit time: 

pUA = m = constant 

(continuity) “ + ~ ~ - 0 (2. 7) 

p u A 

As a result of the assumptions 2 and 3 the equations of momentum and 



state are: 
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(momentum) 




^ = - UdU 
P 


(2.8) 


(state) 


dp 

P 


dp dT dg 

' p T ' (1 - g) 


(2.9) 



From the first law of thermodynamics, with no external heat sources, 
the energy equation is 

d (h + U 2 /2) * 0 . (2. 10) 

The density of the mixture, p, used in the above equations is 
defined by the equation, p = + p'^» The density p'^ is the mass 

of liquid condensate per unit volume of the vapor. This requires an 
assumption that the liquid droplets be uniformly distributed throughout 
the vapor so that the densities of the vapor and of the liquid may be 
referenced to the same volume, g is the mass fraction of the working 
fluid which has condensed into the liquid phase. 

By use of the assumption that the vapor and the liquid are at the 

same temperature, and for constant C and L equation (2. 10) becomes, 

P 

d (U 2 /2 + C T - gL) = 0 (2. 11) 

P 

and 



(energy) UdU + C dT - L dg = 0 

P 



( 2 . 12 ) 



By introducing the sound speed by the perfect gas formula in 

2 R 

terms of vapor density, A = y — T, an alternate form of (2. 12) is 

M 

, n ,,2 dU dT L dg n 

(r-i)M _ + _-_! = o 

P 



(2. 13) 
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The four equations, (2. 7), (2. 8), (2. 9), and (2. 12) involve the six 
unknowns, p, p, T, U, g, and A. For this particular problem, it will 
be assumed that g and A are determined from other sources and p, p, 

T and U are obtained from the above equations. 

Area will be known for a given nozzle, or the area distribution 
must be determined from a given Mach number or pressure distribu- 
tion. The question of determining the condensation, g, depends on the 
assumed flow conditions. If it is assumed that thermodynamic equilib- 
rium exists at all times the Clausius-Clapeyron equation can be used. 
This method will be developed in Section V. However, in the event 
super saturation takes place, the determination of g is not so simple 
and nucleation theory must provide the missing equation. 

After devoting one section to the determination of the onset of 
condensation, and one to nucleation theory, the flow equations will be 
married to the condensation equations for a complete solution. 



IH. PREDICTION OF THE ONSET OF CONDENSATION 



A review of the literature shows that indeed there is no proven 
available theory which is valid for predicting the onset of condensation 
in a supersaturated flow. A solution of the equations for predicting 
this point always requires an assumption as to the amount of condensate 
formed, the droplet formation rate, or the deviation from isentropic 
flow. 

A marriage of the flow equations (thermodynamic) and the conden- 
sation equation from nucleation theory provide a means of computing 
the flow properties for a given area distribution. A plot of a flow param- 
eter during the period of condensation shows a sharp deviation from a 
plot of the same parameter undergoing isentropic expansion. Experi- 
mental investigations have determined that static pressure, p, and pitot 
pressure, p' , are strongly affected by condensation and a jump in these 
quantities usually is used to detect the onset of condensation. 

Here we are defining the point of onset of condensation as that point 
along the flow at which the initial condensate formed is sufficient to 
cause the computed value of p or T to differ a detectable amount from 
its corresponding value for the isentropic expansion. 

The complete solution could be solved for very small steps starting 
at the saturation point and continuing throughout the length of the nozzle. 
From a plot of actual pressure, p, compared to the pressure distribu- 
tion for an isentropic expansion the point of onset of condensation could 
be determined. However, this would require an excessive number cf 
calculations, and for certain combinations of nozzle parameters and 
stagnation conditions the vapor may traverse the nozzle without ever 
condensing and several minutes of computer time would be wasted. 
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Therefore, it would seem expedient to establish a criterion for estimat- 
ing the point of onset of condensation and then to continue the stepwise 
computation of the full set of equations from this point onward. 

Since the expansion of the vapor, in which there is no condensed 
liquid, will be assumed isentropic, any deviation from isentropic con- 
ditions should serve as an indication of condensation. As previously 
stated, the static pressure is sensitive to condensation and this will 
be used as an indicator. From the calculations of Hill (11) the static 
temperature, T, also appears to give a strong indication. At certain 
times it is quite possible that the temperature might give a better indi- 
cation than the pressure, and this factor now will be investigated. 

The flow equations from thermodynamics, (2. 7), (2. 8), (2. 9), and 
(2. 13), relate isentropic conditions when there is no condensate present, 
g = Ag = 0, and relate actual conditions when the vapor is condensing. 

A proper combination of these equations should produce a parameter .. 
which indicates the deviation from the isentrope. Since we are inter- 
ested in temperature and pressure we will recombine the flow equations 




(2. 7) 



dps - pUdU 



( 2 . 8 ) 



dp dp dT dg 




(2. 9) 




T 

P 



(2. 13) 



Substitute equations (2. 7) and (2. 13) into (2. 9) to obtain 



-13- 



dp 

P 



dU 

U 



dA-(v 

A 



1X ..2 dU 
1) M — + 



L dg 

C T 
P 



dg 

(1 - g) 



dp 

P 




dA [ L 

A + |C T 
' P 




J 



(3. 1) 



R 

from (2. 8) and the equation of state, p = p (1 - g) — T, we obtain 



dU 

U 



(1 - g) dp 

P ’ 

y M 



Insert this into equation (4. 1) and rearrange 



dp _ dp j (l - g) 
P P yM 2 
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dA 



C T 1 - g 
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f. 2) 

cr\ 



dg 

j 



For this particular problem we are interested only in the point of 
onset of condensation. Therefore, we assume isentropic conditions 
exist up to the point and the only condensate present will be the dg 
formed over the incremental step at which we are considering the 
onset to take place. The calculations will be accomplished in a step by 
step manner on a digital computer, so we write all increments as A 
and set g = 0 to get 



Ap 

P 



yM 2 
M 2 - 1 



f 




A A L 

+ | cTt 

P 



- 1 



Ag 



(3. 3) 



From (3. 3) it can be seen that over steps in which there is no conden- 
sate (Ag = 0) the stepwise change in pressure will be isentropic. How- 
ever if Ag ^ 0 then the step change in p will differ from isentropic. 

We therefore rewrite (3. 3) as 
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AA 


i A 


L ll 


P 


M 2 - 1 
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Y ~ A A 


C T 

p ' 



Ag, 



(3. 4) 



We now can define 



e = 



p AA 



C T 
P 



Ag 



(3. 5) 



as the fractional deviation from isentropic of the pressure change over 
an increment Ax. The magnitude of which can be considered a sig- 
nificant deviation is yet to be determined. 



Now we examine a similar analysis of the temperature as a possible 
indicator. From equations (2. 13) and (2. 8) we obtain 



= _ g )dp + ^- dg 

T y ' g ' p C T g 



(3. 6) 



P 



By substituting the value of Ap/p from (3. 4) the equation becomes 



AT y - 1 J yM AA 
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Again we define 



e T ' AA 



2 \ 

y - 1/M 



LA 



y - 1 



/ C T 

/ P 



- 1 



Ag 



(3. 8) 



as the fractional deviation from isentropic of the temperature change 
over Ax. 
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Note: A comparison of (3. 5) and (3, 8) shows a difference only of 

the coefficient of L/C^T. Thus as | y | — 1, i. e , as M — 1, 

e m approaches the value of e . However, an examination of e m and e 
T P T p 

for given conditions of L/C T, AA, Ag produces: 

P 



(a) for condensation in the supersonic portion (M > 1) 

(b) for condensation in the subsonic portion (M < 1) — 



_e T >e p 
€ p > £ T- 



As an example we consider reasonable values of L/C T = 3 and M - 10 

A a P 

the result gives = (2) Ag and e T ^ (9. 5) ^Ag, Therefore^ a 

dual criteria will be selected so that will determine the onset of 

condensation when M > 1 and e will be used when M < 1„ 

P 



IV. NUCLEATION THEORY 



A. Condensation Nuclei 

In the process of condensation, the saturated state is referenced 
to a plane liquid phase and the vapor is considered to condense on a 
flat liquid surface or on a "cold" container wall such as the moisture 
formed on the outside surface of a glass of ice water. However, in 
highly expanded flows the liquid surface is not present and the "cold" 
wall is not available, due either to the "hot" boundary layer surrounding 
the flow in a nozzle or to the total absence of a wall in the case of the 
free jet. 

The vapor then must look for other surfaces on which to condense. 
These surfaces appear either in the form of impurities, such as par- 
ticles of dust, etc. , or as small drops formed of clusters of molecules 
of the vapor itself united by statistical collision. 

Obviously, for a truly pure vapor there would be no foreign nuclei; 

3 

however, dust particles in clean air have been estimated upward of 10 
particles per cubic centimeter, so this possibility must be considered. 
Head (2) has shown that these foreign nuclei do play an important role 
in the onset of condensation. However, for highly expanded flows 
where the stay time is very short, Stever (12) has shown that the high 
rate of condensation could not be caused by condensation onto the for- 
eign nuclei alone. Stever (12) quotes an example from Oswatitsch’s 

5 

study of water vapor which considers 10 foreign nuclei per cubic cen- 
timeter traveling in a gas for a distance of 10 centimeters where there 

o 4 

is super cooling of 30 C. A slightly supersonic velocity of 3. 3 x 10 

-4 

centimeters per second would traverse this 10 centimeters in 3 x 10 
seconds. Using the faster growth formula for droplet growth gives a 
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_ 5 

radius of growth of about 3 x 10 centimeters in this distance. The 

-8 

amount of water condensed over this span then is 10 grams per cubic 
centimeter, which releases an amount of heat that is very slight and 
which has an insignificant effect on the flow quantities. The actual 
condensation process has a marked effect upon the flow properties; 
therefore it must be concluded that the number of drops far exceeds 
the number of foreign nuclei. We must then resort to the spontaneous 
formation of the droplets onto which condensate growth occurs. 

There is an opening in the above argument, in that the growth rate 

onto the foreign nuclei could be grossly in error. However, Head (2) 

1 6 

found droplet formation rates on the order of 10 drops per cubic 
centimeter per second in his investigation of water vapor; thus, the 
high particle formation rate tends to support Stever’s argument 

For this analysis it will be assumed that the condensation phenome 
non is one of spontaneously formed nuclei and is triggered by the attain 
ment of critical size droplets onto which the supersaturated vapor 
condenses. 

The concept of the critical drop size as developed in Stever (12) 
stems from the energy required to evaporate a drop of liquid which 
contains n molecules. 

E(n) = nE^ - 4nr 2 o (4.1) 

In (4. 1) is the energy required to remove by evaporation a single 
molecule from an infinite surface of liquid and a is the surface tension. 
To evaporate one molecule would require 

2(J V, . _ 

AE(n) = E - 

00 
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(4.2) 
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where is the volume of a molecule of the liquid. An application 
of Boltzmann's law gives vapor pressure 



P 



K* 



AE 

" kT 



and the equilibrium vapor pressure for a plane surface 
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Combining these gives the Thompson-Gibbs equation 
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Expressing the molecular volume as /VN^p^ and rewriting it as 
a log function, (4. 3) becomes 



In p/p_ = 



2pq 
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oo rp^kN^T rp^RT 



(4. 4) 



Note that with the occurrence of r in the denominator of equation 
(4. 4) the equilibrium between the vapor and the liquid drop is unstable. 
For a vapor which is saturated with respect to a given drcp size, if 
one molecule sticks to the drop the radius increases and the existing 
vapor pressure is greater than the new equilibrium vapor pressure; 
thus, condensation will continue. The reverse is true for evaporation, 
in that if one molecule evaporates the existing pressure is lower than 
equilibrium and evaporation continues. 

Thus arises the concept of the critical drop radius, r% for which 
r < r* the droplet vaporizes and for r > r* the droplet grows. Solving 
(4. 4) gives 
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Hill (11) points out that in a constant environment, if the critical 
radius was used to specify the initial drop size then the growth rate 
would be zero and there would be no condensation because at critical 
size the decay and growth are equally probable. Hill bases his use 
of 1. 3 times the critical radius on Becker and Doering's conclusion 
that at that size the probability of decay is close to zero. In this anal- 
ysis the vapor will be flowing and any critical drops formed will not 
remain in a constant environment very long since the quasi- steady 
steps for the calculations will be very small. Therefore, the com- 
monly used r* of equation (4, 5) will be used as the initial drop size. 



B. Kinetics of Nucleation and Condensation 

As implied previously, the spontaneous nucleation approach will 
be the only one considered in this analysis. There have been volumes 
written in the development of this area and no attempt at reproduction 
will be made here. The general approach and final result of the most 
frequently referenced authors will be stated, however. 



In highly expanded flows the state of the vapor can change very 
rapidly from the unsaturated state to the highly supersaturated state. 
Frenkel (13) treats the equilibrium distribution of drops in an unsatu- 
rated state as a dilute solution of different substances in the vapor as 
a solvent. By classifying the different solutes by the number of mole- 
cules, n, they have in a drop, he obtained the expression for the 



number of drops containing n molecules, 



N = N e 
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where = number of molecules in the vapor state 

^liq' = thermod y nam i c potential of single molecule 

4rrr^a = surface energy term. 

This droplet distribution is based on the assumption of equilibrium 
(<f> y < however, it is generally assumed that the same distribu- 
tion can be used for supersaturated conditions (# < (f > v ). Note that 

for the supersaturated condition the number of drops increases with 
increasing radius and the distribution is unstable. Thus, it becomes 
essential to determine the rate at which the distribution changes from 
the stable distribution of the unsaturated state to the unstable distribu- 
tion of the supersaturated state. 

This requires a determination of the rate of formation of critical 
sized drops, J, in drops per cubic centimeter per second. 



Head (2) states the equilibrium equation of Volmer 

AW 



4tt r^or 



„ kT 3kT 

J = K° e sb K- e 



(4. 7) 



where AW is the total work of the formation of the droplet. 



The determination of the constant, K, has been developed by sev- 
eral authors, including Volmer, Becker and Doring, and Frenkel. 

All assume a quasi- stationary phase transition process in which the 
number of molecules which condense are instantly replaced in the 
vapor such that the number of vapor molecules is maintained constant. 
Volmer considered only the number of drops reaching critical size, 
the number of vapor molecules striking the surface of the critical size 
drop per unit time, and the number of molecules in a drop of critical 



size. 
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Becker and Doring derive the nuclei formation frequency purely 
from kinetic considerations and take into account the evaporation from 
the surface of the drop, as well as the condensation onto it This ap- 
proach is considered by most authors to be better than others; however 
their solution is rather difficult to apply. 

Frenkel uses the method of Zeldovich to integrate the differential 

equations of Becker and Doring and gets a much simpler derivation of 

an equation that is almost identical to that of Becker and Doring. The 

2 3 

ratio of the result of Frenkel to that of Becker-Doring is (n*) y and 
this is considered minor, n* is the number of molecules in a critical 
sized drop. 

A partial derivation of the nucleation rate equation is presented 
in Appendix A. This derivation is necessitated by an apparent error 
of 7T in the original equation presented by Frenkel (13). 

For this analysis, equation (A. 7) will be applied as the nucleation 
rate equation. 



where J = number of critical drops formed per unit volume per unit 
time. 

C. Limitations 

1. Generally a as a function of temperature is not known and this 
is a large factor in the result. Surface tension appears to 
the third power in the exponent; hence a small error will have 
a large effect upon J. 
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2. The rapidity of expansion will show marked deviations from 
the assumed quasi- stationary condition. 

3. The condensing flow may have sublimation directly from vapor 
to solid which is not considered. 

4. At high degrees of super saturation the drop radius will be small, 
and Head (2) has shown that the variation of a with r is large 

for drops of few molecules. Stever and Rathbun (5) have included 
this effect into the Frenkel equation to obtain 



Bogdonoff and Lees (8) have an approach which is entirely dif- 
ferent from that of Head and Stever and which produces quite 
different results. Until these theories have been developed 
further, it is not practical to include either theorem in this 
analysis. 



From the previous discussion of the assumptions made in the formu- 
lation of the nucleation theory and of the limitations above, it is apparent 
that the best one can expect in an analysis such as this will be an order 
of magnitude estimate of the condensation parameters and a qualitative 
analysis of the flow conditions. 

D. Growth Estimates 

Heretofore, in order to calculate the formation of the critical drops, 
we have considered the temperature of the liquid drops to be the same 
as the surrounding vapor. In actuality, the impinging and vaporizing of 
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molecules onto and away from the droplet leave a surface temperature 
that approaches the saturated equilibrium value for the surrounding 
vapor. This temperature difference allows for growth of the particle. 
The growth of small drops, i. e. , r << mean free path, occurs by the 
molecular transport of molecules onto the surface, some of which con- 
dense, and by the conduction away of latent heat of vaporization by the 
evaporation of a fewer number of molecules which leave at higher kinetic 
energies. 

Oswatitsch laid the basic groundwork for growth in a flowing vapor. 
Stever (12) follows the method of Oswatitsch to develop aq equation for 
the growth rate of the small drop. 



The development presented by Stever makes use of the kinetic 
theory of free molecule flow (See Patterson (15)). Consider a drop 
of radius r, temperature T , and saturation pressure (ps^)^ with 
a surrounding vapor of pressure p and temperature T. The heat trans- 
fer from drop to vapor per unit time is 



P 



Sff (T liq ' T > “ 



where a is an accommodation coefficient which varies from 0 to 1, 
depending on the vapor and liquid surface. 



Also from kinetic theory the heat transferred to the drop surface 
due to condensation is Lfp 7 kT/ 27 r m where p 7 kT / 2 tt m is the amount 
of vapor (in the form of molecules) impinging on a unit surface and f is 
the fraction which condenses. The energy being transferred to the sur- 
face is equated to the energy being transported away to obtain an expres- 
sion for f; 
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the mass of the drop is then increased by 
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If the calculation is made in quasi-steady steps, time can be expressed 
as a function of step length and 
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(4. 10) 



where the accommodation coefficient, a, is a constant, approximately 
equal to 1 and is the temperature of the surface of the liquid drop. 

The accommodation coefficient a is expressed by Kennard (16) as 

a (E. - E ) = (E. - E ) 
l w l r 

where E^ - Energy brought up to a unit area per unit time by the 
impinging molecules 

E^ - Energy carried away from the molecules leaving a 
unit surface per unit time 

E^ = Energy that the departing molecules would carry away 
if they carried away the same mean energy per mole- 
cule as does a stream issuing from a gas in equilibrium 
at the surface temperature. 
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Patterson (15) lists coefficients of the range from . 8 to 1. 0, 
and Stalder (17) found that for free molecule flow an accommodation 
coefficient of 0. 9 checks well with test results. Thus we arrive at 
the statement, "approximately equal to 1. " This only applies for cases 
where free molecule flow can be assumed to prevail. 



On the other hand, where r approaches the magnitude of the mean 
free path the heat removal is controlled by normal gaseous heat con- 
duction. For this large drop Stever (12) gives the growth equation for 
period of time, t, 
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Other authors have made similar analyses, and Wegener (14) re- 
cords the quasi-steady solution of Buhler (18) which includes an expres- 
sion for the condensation coefficient. 
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(4. 12) 
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where T g , the surface temperature of the drop, may be chosen as the 
equilibrium saturation temperature of the vapor with respect to the 
drop. As in the previous case, Wegener states (4. 12) is valid only 
for small drops and the equation recommended for large drops is 
merely a revised form of (4. 11). 



Hill (11) notes that typical calculations of condensation nuclei indi- 
cate very small sizes which are generally much smaller than the mean 
free paths; thus only the free molecule equations need be considered. 
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His analysis produces a growth rate expression of 



f = [ rt d - RT 1 

2tT, p T / 2 tt RT 
fg K L 



(4. 13) 



in which - change of internal energy due to condensation. The 
term is droplet temperature, and Hill presents calculations to show 
that the rate of drop temperature change is enormous compared to the 
rate of growth. For water vapor the temperature approaches wet bulb 
temperature in approximately the same time that it takes the number 
of molecules in the drop to increase by a factor of 10. Hill uses this 
in expaqation of the assumption used by Oswatitsch that the drops al- 
ways have wet bulb temperature during the growth period. 



After a displacement, the drop temperature will rapidly approach 
the temperature at which d T^dt - 0 at a much greater rate of change 
than dr/dt . Hill uses this to form an expression for which depends 
only upon the instantaneous values of r, p and T. 
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(4. 14) 



where £ is the fraction of impinging molecules which condense. 



From the investigations on water vapor, nitrogen, air and a few 
metal vapors it appears that for highly expanded flows the small drop 
size is the general rule. The calculations will then be initiated by 
applying (4. 10) as the growth rate equation and choosing drop tempera- 
ture as the equilibrium saturation temperature of the vapor. The ac- 
commodation coefficient in (4. 10) will be varied to study the effect and 
a will be placed into the computer program as a data input so that accu- 
rate values, when known, can be used for the calculation. 



V. ANALYSIS 



The basic equations are established in Sections II through IV. The 
proper combination of these equations produces a theoretical solution 
for the hypersonic expansion of a pure vapor in which a supersaturated 
state occurs. Another solution is obtained for which there is saturated 
equilibrium expansion in lieu of the supersaturated portion. Both solu- 
tions, however, stem from the same basic model. 

The basic model considered here consists of assuming: 

1. "One-dimensional nozzle" with geometrically specified inlet 
and diffuser. 

2. A pure vapor with specified stagnation conditions and known: 
constant specific heats, latent heat of vaporization; surface 
tension, for infinite radius as a function of temperature; and 
saturation curve. 

3. An isentropic expansion to the point of onset of condensation 
and then diabatic flow through the region in which condensa- 
tion is taking place. 

A. Basic Assumptions 

In addition to the assumptions of the basic model, this analysis 
also embodies the multitude of simplifying assumptions which have 
been included in the derivation of the basic equations. The type of 
equations and the specific assumptions made are: 

1. Flow equations 

a. Mass flow is constant 

b. Flow is one -dimensional and steady state 

c. Volume of condensed phase is negligible when compared 
to the total volume 
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d. Nozzle is frictionless with no heat transfer across nozzle 
wall. 

e. Condensed mass is liquid, is uniformly distributed through- 
out the gaseous components, and has same speed and tem- 
perature as the stream 

2. Prediction of onset of condensation equations 

a. This embodies same assumptions as flow equations 

b. No condensate is formed prior to the predicted point of 
onset of condensation 

3. Nucleation equations 

a. Drops are formed by spontaneous nucleation only 

b. Only drops reaching critical drop size continue to grow 

c. Number of molecules of vapor is maintained constant 

d. Saturation pressure corresponds to saturation pressure 
of a droplet of infinite radius 

e„ The ordered velocity of the particles is neglected 

f. Equilibrium values of the probabilities that one molecule 
will leave from or condense on a unit surface and equilib- 
rium particle distributions are assumed valid for the 
non -equilibrium case 

4. Growth rate equations 

a. Drop radius is much less than mean free path 

b. Accommodation coefficient, a, is approximately equal to 

1 . 0 

c. Drop temperature is assumed equal to the saturation tem- 
perature that corresponds to the pressure of the vapor 
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B. Supersaturated Expansion 

The supersaturated expansion passes the saturation point without 
any condensate being formed; therefore, a four step solution is required: 
1. Isentropic equations from Shapiro (19) and NACA Report 1135, 
are used to calculate the expansion of the vapor to the point of 
condensation. 
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( 5 . 1 ) 



(5. 2) 



(5.3) 



(5. 4) 
(5. 5) 



At any particular value of T the above equations are readily 
evaluated for given p q , T q , p q , 7, and A*. The station at 
which these values apply is then determined from the given noz- 
zle geometry^ 

2. The saturation point is established by an iterative balance of 
the saturation curve equation with the isentropic flow equation. 
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The saturation curve can generally be expressed in the 

c 

form In p = B - — , in which B and C are constants whose 
oo T 

value depend upon the choipe of units for p and T. The equa^ 
tions which require joint solution are 
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The isentropic relations above are then utilized to calcu- 
late M P ,, T ., A , . p ,, and U ,. 
sat sat sat’ sat, sar sat 

3. The point of onset of condensation is determined by taking 
steps, AT, of super saturation and computing the flow param- 
eters from the isentropic relations for each successive step. 
The isentropic values of each incremental step in temperature 
are introduced into the nucleation equations for critical drop 
size and formation rate. 
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(4. 8) 



Note that choice of a very small AT will cause the first 
step to remain very near the saturation point and the resultant 
p/p^ = 1 will cause r* to approach an infinite radius, r* falls 
off very rapidly, however, as the degree of super saturation 
increases, and for most cases a choice of AT > 5°k is sufficient. 
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The growth rate is zero for the critical droplet at the instant 
of formation; therefore the nucleation rate for this increment is 



N T ' J T A T AX T 



(5. 8) 



The mass fraction of condensate formed in this increment is 
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(5.9) 



The epsilon equations from the onset of condensation develop- 
ment, (3. 5) and (3. 8), are now used to determine if this amount 
of condensate is sufficient to cause an appreciable deviation 
from the isentropic expansion. 
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If the e is less than the assumed critical value then the same 
calculations are carried out at step T + AT until the critical 
value of e is equalled or exceeded. If exceeded, a bracketing 
procedure is employed to improve the estimate. 

All quantities computed for Ag at values of e < € cr ^ ra j 

are discarded, and it is assumed that there has been no con- 

densate formed prior to e - e ... 

critical 

The choice of Ax in equation (5. 8) will effect directly the 
magnitude of N and of Ag in equation (5. 9). AA in equations 
(3. 5) and (3. 8) is also a function of Ax so that in the evalua- 
tion of e the effect of Ax is nullified. Calculations have been 
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performed for values of Ax ranging from . 01 to 1. 0 cm to verify 
that the quantity has no effect on the prediction of the point of 
onset of condensation. Since the Ax that will be used for incre- 
menting the calculations through the condensing portion of the 
nozzle will be a required input to the computer program, it is 
advantageous to use that same Ax in evaluating (5. 8). 

4. The condensing portion of the flow requires a joint solution of 
the nucleation and growth equations and the diabatic flow equa- 
tions. The calculations are performed for increments of Ax 
starting from the onset of condensation. Values are estimated 
of p, T and U at the point i and Ag^ is calculated from the 
equations 
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The subscripts i and j are numbered increments of x 
starting from the condensation point, i acts as a label for 
each group of droplets of a particular size and denotes the 
particular increment in which these drops originated as crit- 
ical sized drops, j denotes the increments presently under 
consideration. Thus a designation denotes droplets 
which were formed as critical drops in increment 1, have 
undergone the growth Arg, and are now undergoing the 
growth Ar^. 

Using the above approximate values of Ag and g and the 
unknown values of AA and A the equations from Section n 
are solved for p, T, p and U. 
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The solution of the flow equations is accomplished by the 
following iterative procedure: 



a. 

b. 

c. 

d. 

e. 



+ AU 



assume AU and compute U. - U. 

J J " 

solve (2. 12) for AT and compute T^ - T^ 
solve (2. 7) for A p/p 



AT 



solve (2. 9) for Ap/p 

solve (2. 8a) for AU and check with value from step (a). 
Improve the assumed AU and continue until a solution is 
obtained. 
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These corrected values of p, T and U are now used to 
resolve the equations for Ag, and the process is continued 
until a set of flow parameters which satisfy both sets of 
equations is obtained. 

This iterative procedure is continued for each increment 
Ax up to a specified x max > usually the end of the nozzle. 

This set of equations could be evaluated from the satura- 
tion point onward in order to obtain a more accurate predic- 
tion of the onset of condensation. It is obvious, however, that 
the magnitude of the calculations involved is formidable; and, 
the slight error induced by attempting to predict the point of 
onset of condensation, does not justify the refinement. This 
is especially true since the approximations in the basic 
equations are so gross. 



C. Saturated Equilibrium Expansion 



The saturated equilibrium expansion is a gradual process with 
condensation starting at the saturation point and equilibrium flow 
being maintained throughout. This occasionally occurs for slow 
expansions of certain vapors and for cases in which the impurity 
content is high. In this analysis the saturated equilibrium expansion 
is presented only as a comparison with the supersaturated case. 

The calculations are greatly simplified in that the mass fraction 
of condensate can now be specified by the Clausius -Clapeyr on 
equation 
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This equation and the flow equations (2. 7), (2. 8), (2. 9) and (2. 12) 
usually are integrated and expressions are written for each param- 
eter (see Wegener (14) or Daum (10)). By proper manipulation and 
combination the integrated equations can be reduced to one involving 
g. 



g = 






(5. 13) 



However, in this analysis the flow equations are programmed in the 
incremental form for the supersaturated case, so it is convenient to 
write (5. 13) in the form 



Ag 
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(5. 14) 



and to use this for the Ag expression in a solution similar to that 
used for super saturation. 

1. The saturation point is established in the same manner as 
for the supersaturated case. Small increments of length, 

A x, are started at the saturated point. 

2. An incremental increase in the mass fraction of condensate, 
Ag, is assumed for the step Ax. With assumed value of Ag 
and known value of A the flow equations are solved for p, p, 
and T. 
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Equation (5. 14) is solved to check and correct the assumed value of 
Ag and an iterative procedure is established to obtain a satisfactory 
solution of the five equations. Again this is continued in a step by 
step process until the solution is complete., 
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VI. DIGITAL COMPUTER PROGRAM 



A major task in this study has been the adaptation of the mathematical 
analysis of Section V to the digital computer. The description contained 
herein is presented for the future user of the program and in no way 
covers the myriad of detours encountered between the input and output 
statements. 

The Computer Center of The University of Michigan has developed 
a Michigan Algorithmic Decoder (MAD) computer program which com- 
piles much faster than the commonly used FORTRAN program. The 
computer program for this investigation has been written in the MAD 
language and can be translated into FORTRAN if desired. The IBM 7090 
computer program for the condensation of a pure vapor is presented in 
Appendix B. 

The program has been subdivided into a main program and six 
subroutines. The subroutines enable easier modification of the program 
and provide a greater flexibility in its use. Statements as to the func- 
tion of each of the calculations, the input parameters, and the obtainable 
results are provided in the program at the beginning of each subroutine. 
Appendix B also includes a list of the program symbolic names assigned 
to the parameters used in the equations of Section V. 

A brief summary of the subroutines follows: 

1. Isentropic flows (IFLOWS) — Given T, this program solves 
the isentropic flow equations for p, p, U, M, and A. 

2. Condensation Flow (CFLOW) — For given g and A, this pro- 
gram solves the diabatic flow equations (2. 7 through 2. 12) 
for p, T, p, M, and U. 
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3. Nucleation Theory I (NUCLEI) — From nucleation theory equa- 
tions (4. 5, 4. 8, 5. 8 and 5. 9), NUCLEI computes the condensate 
formed due to the nucleation of critical sized drops only, i. e. , 
growth is not included in this portion. This subroutine is used 
to determine condensate formed for that portion of the main 
program which determines the point of onset of condensation. 

4. Nucleation Theory II (NUCLE2) — -This subroutine evaluates 
the growth equation (4. 10). NUCLE2 is used in conjunction 
with NUCLEI to determine the amount of condensate present 
for all calculations downstream from the point of onset of 
condensation. 

5. Nozzle (NOZZLE) — The geometry of the nozzle is placed in 
this subroutine, from which A, AA, or x can be determined. 

6. Vapor Parameters (FRHOL, FL, FSIGMA) — The parameters 
p^, L and a are placed in this subroutine as functions of T 

or as constants and are evaluated as called upon by the 
Main Program. 

The Main Program consists primarily of a proper utilization of 
the subroutines. The operations by the Main Program appear in the 
following order: 

1. Read in the input data. 

2. Use IFLOWS to expand isentropically from stagnation condi- 
tions. 

3. Match isentropic equations to saturation curve to determine 
saturation point. 

4. Use IFLOWS for expansion past the saturation point and 
calculate Ag from NUCLEI. 

5. Evaluate epsilon equations (3. 5 or 3. 8) to determine onset 
of condensation. 
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6. Downstream of the onset of condensation, a joint solution of 
NUCLEI, NUCLE2, and CFLOW is obtained at each incre- 
ment, Ax, and the parameters x, N, Ax, A, AA, P, p, T, 

U, M, r*, r.„, Ar, Ag and g are stored for later print out. 

The saturated equilibrium expansion is compiled in a very simple 
program and need not be presented at this time. 

The Nozzle and Vapor Parameter Subroutines are included in a 
Program Common such that any of the included parameters, when used 
in the Main Program, are automatically evaluated and returned. 

Conditional statements are inserted in the program to account 
for the possibility of condensation occurring in either the converging 
or diverging portions of the nozzle, the possibility of reaching the end 
of the nozzle without any condensation having occurred, and the failure 
of the iteration of NUCLEI and NUCLE2 with CFLOW to converge. 

The development of the equations is such that any consistent set 
of units can be used for the problem. However, it was found that 
minimum scaling was necessary if the problem was worked in a cgs 
system of units. The problem, as written, is compatible with the 
cgs system of units and a list of the proper input data is presented 
in Appendix C. 



Vn. CALCULATIONS AND RESULTS 



The calculations performed in this study are divided into four 
parts: 

1. A test of the program by hand calculations and a test of the 
method by comparing the results for a nitrogen expansion 
with the experimental results of Willmarth (3). 

2. Justification of the assumptions of small radii, no drag, and 
accommodation coefficient approximately equal to one in the 
basic assumptions. Verification of the prediction procedure 
for the onset of condensation. 

3. A variation of latent heat, specific heat, surface tension, and 
nozzle geometry for the nitrogen expansion to observe the 
induced effect of each parameter. 

4. Application of the theory to metal vapors by computing the 
expansion from an initial pressure of approximately one 
atmosphere of copper, tin, lead, and zinc vapors. 

A. Validation of the Program 

The input data for the calculations performed in this study are 
tabulated in Table C-l of Appendix C. A reference to Nitrogen (1), 
Zinc (I), etc. , indicates the set of input data from which the indicated 
results were obtained. Any additional modifications of the input data, 
such as a change in a, L, or nozzle geometry, are listed directly on 
the plot of results. An explanation of the input data required, the 
proper units for each parameter, and the references (20-24) from 
which the vapor properties were obtained is also presented in 
Appendix C. 
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Initially the computer program was divided into small segments 
and hand calculations were performed for Nitrogen (I) to verify the 
computer results. With the assurance that each small segment was 
working properly, the entire program was assembled and hand calcu- 
lations were performed to spot check the calculations of the saturation 
point, the point of onset of condensation, and the first, few incremental 
steps of the condensing flow regime. 

With this degree of confidence in the program, the next step was 
to check the theoretical results with known experimental results. The 
Nitrogen (I) input data is designed to approximate an experimental 
test of the condensation of nitrogen performed by Willmarth (3). 

Nitrogen properties from NBS Circular 564 (24) and the Interna- 
tional Critical Tables (23) required extrapolation below 70°Kfor the 
surface tension and the latent heat of vaporization. These values are 
generally valid only to within ±10% and, as will be shown later, this 
becomes extremely critical in the case of the surface tension. The 
vapor pressure curve for nitrogen is reasonably well known, and it 
was approximated by the Clausius-Clapeyron equation (see Fig. C-l). 

For comparison with experiment the calculations for nitrogen were 
carried out in full from the saturation point onward and no attempt was 
made to predict the point of onset of condensation. A comparison of 
the experimental results of Willmarth (3) to the theoretical prediction 
of the condensation of Nitrogen (I) is presented in Figure 1. The the- 
oretical curve for Nitrogen (I) shows a degree of supersaturation of 
approximately 3. 5 degrees (20%) lower than is indicated by experimental 
results. This is expected in that the amount of impurities in the experi- 
mental nitrogen was unknown and this advanced nucleation would cause 
an early break. However, of greater importance is the degree of 
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uncertainty in the extrapolation of the available surface tension data. 

In anticipation of the effect of a on the theoretical results, a second 
extrapolation was chosen for a. The a equation for Nitrogen (II) re- 
mains within the ± 10% accuracy of the original surface tension data 
and, in actuality, gives a value for a which is only 7% less than the 
ct from Nitrogen (I) at the onset of condensation. 

The theoretical prediction for the condensation of Nitrogen (II) 
compares most favorably with the experimental results in Figure 1. 
The strong influence shown by only a 7% change in a will be discussed 
more fully at a later point in the discussion. Figure 2 presents a 
more complete comparison of the theoretical results for Nitrogen (II) 
as compared to the experimental results of Willmarth (3). The isen- 
tropic and saturated equilibrium expansions are also presented in this 
figure. The theoretically calculated results presented in Figure 2 
employ the full calculation of the nucleation and flow equations for the 
entire length of the nozzle. This is significant in that the flow prop- 
erties prior to the "condensation shock" lie on the isentrope and the 
flow properties after the "shock" tend to approach the saturated 
equilibrium expansion. This is the expected result, and it gives 
further confidence in the program. 

The experimental points of Willmarth (3) fall on the isentrope as 
well as along the "condensation shock", which indicates a proper esti- 
mate of the physical expansion and validates the program. The close 
correlation between theoretical and experimental results within the 
"shock” itself lends further support to the basic nucleation theory of 



the condensation process. 
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Willmarth (3) obtained the pressure results, p, p q and p Q f by 
direct measurement. The effective area ratio, (A/A*)^, tempera- 
ture, T, and mass fraction of condensate, g, were then calculated 
from the measured quantities by a joint solution of the flow equations. 

B. Justification of Assumptions 
1. Drop Size 

A typical distribution of maximum and minimum droplet 

radius is presented in Figure 3 for the Nitrogen (I) expansion. 

Note the very short span over which the nucleation of critical 

sized drops, J curve, takes place. The average critical sized 

-8 

drop has a radius of 4. 5 x 10 cm, and the largest particles 

-7 

which leave the nozzle have a radius of 8 x 10 cm. The radius 

of the nitrogen molecule from Loeb (25) is approximately 
-8 

1. 5 x 10 cm and calculated mean free paths vary from the 

-4 -3 

order of 10 cm at the onset of condensation to 10 cm near 

the nozzle exit. According to Stever (5) the drop radius is 

larger than the minimum required to contain the 10-12 molecules 

necessary for the application of macroscopic theory. 

Therefore it is not unreasonable to assume that small drops 
are the general rule, that the no drag approximation is valid, 
and that since r << mean free path the samll drop growth equa- 
tion can be used. 

2. Accommodation Coefficient 

The effect of a variation in a is presented in Figure 4. Small 
variations near the value of a = 1. 0 have very little effect. How- 
ever, since the growth equation is a direct function of a, the curve 
for a - 0. 1 differs markedly from a = 1. 0. Fortunately, the known 
values for a are from 0. 8 to 1. 0 so that the choice of a = 1. 0 is a 
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reasonable one. The computer program is so constructed that 
accurate values of a, when known, can be inserted as input 
data. 

An additional feature of Figure 4 is the indication that growth 
has little effect up to the "break" point but thereafter controls the 
remainder of the condensation process. This is further brought 
out by the plot of J in Figure 3, which indicates the rapid nuclear 
tion of a large number of particles to start the "shock" and then 
a continuing growth of these particles until saturated equilibrium 
conditions are met. 

3. Point of Onset of Condensation 

From the previous discussion and from Figure 3, it appears 

that little or nothing happens in regard to condensation until a 

degree of super saturation is reached for which the critical drop 

radius is very small. As r* from equation (4. 5) becomes small 

2 

the exponential term, which contains r* , ceases to be the pre- 
dominant term in equation (4. 8) and J increases rapidly. 



, 2j U£ 

1 p RT In p/p 
"L oo 



(4. 5) 



J = 




. *2 
4tt ar* 




(4.8) 



The epsilon of Section III should then predict the approximate 
point of onset of condensation. Several calculations for various 
values of € are presented in Figure 5. e = 0 indicates that the full 
calculation of the nucleation and flow equations was carried out 
from the saturation point onward. The choice of epsilon results 
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in only a minor error when compared to the possible 10% error 

-3 -4 

in surface tension. Epsilon values of 10 and 10 give almost 
identical results, thus e = . 001 appears to be a suitable choice 
for the expansion of nitrogen. 

Also in agreement are Head’s (2) calculations for water vapor 

1 6 

in which he indicates that the attainment of J = 10 droplets per 
centimeter cubed per second signifies the point at which the con- 
densation shock will break away from the isentrope. Several 

1 6 

calculations on the computer indicate that J = 10 at the break- 
away point for nitrogen also. At present, however, this author 
sees no justification for predicting that this will be true for other 
vapors as well. 

C. Variation of Parameters 
1. Surface Tension 

As indicated previously, the variation of a is the most critical 
aspect of the theoretical prediction of the condensation phenomenon. 
From equations (4. 5) and (4. 8) it can be seen that a enters the 
exponential term of J as a third power. Thus for a 10% change 
in a the nucleation rate changes by several orders of magnitude. 
Figure 6 and Figure 7 present several variations in a to show the 
large effect that this parameter has on the degree of super satura- 
tion. Later results will show that the effect of other parameters 
is minor when compred to surface tension. 

There appears here a possible method of improving the surface 
tension data. The results of a closely controlled experimental vapor 
expansion could be used as a reference plot and surface tension 
varied on the computer until a matching plot is obtained. Providing 
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the other parameters are reasonably well known, this should pro- 
vide a better estimate of a than merely extrapolating a curve far 
away from the few experimentally determined values that are 
available. 



2. Specific Heat 

0^ is reasonably well known and only varies a few percent 

over the range considered in this analysis. For this analysis 

C is chosen as a constant, and from the plot in Figure 8 the 
P 

constant assumption for appears reasonable. The effect of 

a 10% variation in C is negligible when compared to an equiva- 

P 

lent uncertainty in a. 



3. Latent Heat of Vaporization 

Latent heat is generally a weak function of temperature. For 
Nitrogen (I) the variation in L is only 9% over the range of tem- 
perature from saturation to the nozzle exit. From the plot of 
variations in L presented in Figure 9 it is evident that the con- 
stant latent heat assumption in the basic equations is a reasonable 
assumption and any variation in L is negligible when compared 
to the effect of an equivalent variation in a. 



4. Nozzle Geometry 

The condensation process is a function of time, so it is ex- 
pected that the rate of expansion should have some effect on the 
’’condensation shock. ” Theoretical calculations were performed 
for the Nitrogen (I) expansion in a two dimensional nozzle for which 
the exit half-angle was varied from 9. 75 to 45 degrees. Nozzle 
lengths were selected so that a Mach 8 nozzle was used for each 
run. The results of varying the rate of expansion for Nitrogen (I) 
are presented in Figure 10 and Figure 11. In Figure 10 it is noted 
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that the faster expansion causes the vapor to supersaturate to a 
lower temperature before the onset of condensation. The "conden- 
sation shock" is thinner for the fast expansion, it is (in physical 
dimensions) 3 cm for the 9 - 45° expansion versus 5 cm for 

the 9 ^ = 9. 75° expansion. 9 is the nozzle exit half -angle. 

This is somewhat misleading because, for Figure 11, in which 
all nozzles are referenced to a non-dimensional length for a Mach 
8 nozzle, the faster expansion tends to smear out over a broader 
relative length. The plot of the percent of condensed vapor in 
Figure 11 indicates that the higher the rate of expansion the less 
moisture is formed in the nozzle. It appears that for sufficiently 
high rates of expansion a condensation free flow can be obtained. 
This has been borne out by "condensation free" experiments in 
the past. 



D. Application to Metal Vapors 

As stated in Section V, the theory as developed herein should apply 
to metal vapors. The primary interest here is in the finding of a metal 
vapor that will condense in a highly expanded flow and that is suitable for 
laboratory experimental use. Hill (11) has predicted that sodium and 
potassium will condense; however, these vapors introduce critical han- 
dling problems. Copper, zinc, tin and lead appear to have possible experi- 
mental use, with copper being the preferred vapor of the four. 

Theoretical calculations were performed for a two-dimensional nozzle 
with a constant slope diffuser for the metal vapor input data listed in 
Appendix C. 
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1. Copper 

Fortunately, copper vapor is theoretically predicted to condense 

for a fairly rapid expansion in a two-dimensional nozzle with a 20 

degree diffuser half-angle. The results are presented in Figure 12 

for an expansion in a Mach 10 nozzle from an initial pressure of 

approximately one atmosphere and an initial temperature of 3000°K. 

The "condensation shock" occurs very near the throat, and the change 

is very rapid with the "shock" spanning less than a centimeter of 

nozzle length. The average critical drop radius is approximately 

3 x 10 cm and the maximum droplet formation rate is of the order 
19 

of 10 drops per centimeter cubed per second. Again the onset of 

1 6 

the condensation appears to be marked by a value of J = 10 or 
17 

10 drops per centimeter cubed per second. Copper could be a 
suitable choice for experimental work in rapid expansions because 
the large fraction of condensate formed in this theoretical solution 
indicates that higher rates of expansion and higher initial tempera- 
tures and pressures probably could be used without resulting in a 
"condensation free" expansion. 

2. Z inc 

Two computer runs were made for zinc, and neither run pre- 
dicted any detectable amount of condensate. Zinc (I) was computed 
for a 45 degree diffuser half-angle and initial pressure and tempera- 
ture of approximately one atmosphere and 1500°K respectively. The 
results of this calculation are presented in Figure 13 where a suffi- 
cient number of the calculated points are plotted along the isentrope 
to show that the expansion is essentially isentropic. This indicates 

a negligible amount of condensation. The tabulated calculations from 

~8 

the computer give a maximum g of 1. 48 x 10 at the nozzle exit. 



-49- 



The maximum nucleation rate obtained is J = 2 x 10 drops per 

centimeter cubed per second and the critical drop radius is 
-8 

2x10 cm. It was anticipated that a slower expansion with a 
lower starting temperature would predict a detectable amount of 
condensate. Zinc (II) was run on the computer for a 20 degree 
diffuser and a pressure and temperature of one atmosphere and 
1300°K respectively. The results in Figure 13 again plot along 
the isentrope and indicate no condensate. The maximum g of 

_5 

4. 8 x 10 for the Zinc (II) run indicates some movement toward 

14 

condensation. The maximum nucleation rate of J = 3 x 10 drops 
per cubic centimeter per second is obtained at a vapor temperature 

Q _ Q 

of 375 K and the critical drop radius is again 2 x 10 cm. The 
14 

J - 10 value with no condensate formed lends support to the 
1 6 

J ~ 10 criterion for the onset of condensation. A slower expan- 
sion of zinc probably would produce a condensation shock; however, 
this would be out of the desired rapid expansion regime. From 
the results of the Zinc (I) and Zinc (II) runs it appears that for 
rapidly expanding flows zinc vapor may not condense. 

3. Tin and Lead 

Unfortunately, the computer runs for tin and lead failed to 
converge in the CFLOWS Subroutine. The reason for this is not 
immediately apparent, and a further investigation will have to be 
made before any definite conclusions can be drawn. However, the 
few points that were calculated show promise that both tin and leaf! 
will condense in a rapid expansion. 



vm. CONCLUSIONS 



An adaptation of the digital computer to the analysis of the condensation 
of a flowing vapor is possible, and the program for a pure vapor is pre- 
sented herein. Theoretical predictions of the condensation of nitrogen, 
copper and zinc indicate the following: 

1. The nucleation theory of Frenkel (13) provides a reasonable 
approximation to the rate of formation of condensate in a flow- 
ing vapor. 

2. The theoretical prediction for the condensation of nitrogen 
compares surprisingly well with the experimental results of 
Willmarth (3). 

3. A small drop size is the general rule. The critical drop radius 

-8 

is of the order of 10 cm and droplets of maximum growth 
“ 6 

rarely exceed 10 cm. 

4. The nucleation rate controls the onset of condensation, and a 
prediction procedure based on this parameter might prove more 
useful than the one proposed herein. 

5. The growth equation plays a predominant role once condensation 
has started and a more accurate expression would be beneficial. 

6. The surface tension is the most influential parameter. Generally 
the value of o is uncertain and this leads to the largest source 

of error in the program. 

7. Values of specific heat and latent heat of vaporization are not 
critical, and treating these parameters as constants is a rea- 
sonable assumption. 

8. The condensation process is very sensitive to rates of expan- 
sion, and for rapidly expanding vapors a "condensation free" 
expansion is possible. 
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9. Copper vapor will condense in rapidly expanding flows and 
appears suitable for experimental work in the condensation 
of metal vapors. 

10. Zinc vapor does not condense for high rates of expansion from 
one atmosphere pressure, and runs for lower pressures and 
temperatures are recommended. 
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Figure 1. Theoretical Results for Nitrogen Compared 
to Experimental Results of Willmarth (3). 
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Figure 2. Theoretical Prediction of the Condensation 
of Nitrogen (II) Compared to Experimental Results. 
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Figure 3. Typical Expansion of Nitrogen (I) Presenting Variation of 
Nucleation Rate, Maximum Droplet Radius and Critical Droplet Radius 

Along Nozzle Axis. 



Drop Radius, 



-58- 




b£ 




x, cm 

Figure 4. Theoretical Prediction of the Effect of a Variation in the 
Accommodation Coefficient on the Condensation of Nitrogen (I). 
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Figure 5. Effect of Choice of e as a Predictor of the Onset of 
Condensation for Nitrogen (T). 
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Figure 6. Effect of ± 10% Variation in Surface Tension on the Theroetical 
Prediction of the Condensation of Nitrogen. 
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Figure 7. Effect of Variations in Surface Tension on the Predicted 
Supersaturation of Nitrogen. 
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Figure 8. Comparison of 10% vid 50% Changes in Specific Heat to a 10% 
Change in Surface Tension for Nitrogen (I). 
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Figure 9. Comparison of 10% ana 50% Changes in Latent Heat of Vaporization 
to a 10% Change in Surface Tension for Nitrogen. 
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Figure 10. Theoretical Expansion of Nitrogen (I) for a Two-Dimensional 
Nozzle with Varying Exit Half -Angle, 9 
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Figure 11. Theoretical Expansion of Nitrogen (1) in Mach = 
8 Nozzle of Lengths 35, 16. 5, and 6 cm. 
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Figure 12. Theoretical Prediction of Condensation for the Expansion 
of Copper Vapor in a Two-Dimensional Nozzle. 
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Figure 13. Theoretical Prediction of the Expansion of Zinc (I) 
and Zinc (II) Vapor Indicating No Condensation. 



Appendix A 

DERIVATION OF THE NUCLEATION RATE EQUATION 



This is in no way intended to be a complete derivation of the nuclea- 
tion rate equation. The full derivation can be found in complete detail 
on pages 366 through 400 of Frenkel (13). A scant derivation is presented 
in Stever (12), pages 532 to 540. 

There appears to be an error of n in the final equation presented 
by Frenkel, so enough of the derivation will be reproduced here to clear 
up this discrepancy. The solution is a statistical problem of kinetic 
theory. 



In an equilibrium distribution of droplets where 
lowing relation must hold. 



<</>.. the fol- 
hq 



N S a = N , S 
nnn n-ln-1 




(A. 1) 



where <t>^ and are the thermodynamic potentials of the vapor and 
liquid drops, respectively. N^ and N^ ^ are the number of droplets 
containing n and (n - 1) molecules, respectively. and are 

the surface areas of droplets containing n and (n - 1) molecules, respec- 
tively. is the probability that one molecule leaves a unit surface 
area of a droplet containing n molecules. 



B „ is the probability that one molecule condenses on a unit surface 
n - 1 

area. Assume a and B -.to remain valid for the non-equilibrium case 
n n - 1 

(<£ < (j > v ) and modify the distribution N^ to a non-equilibrium distribu- 

tion f . The equation will no longer balance so (A. 1) becomes 
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J = f , S 
n n - 1 n 



1 - 1 



- f S a ^ 0 
n n n 



(A. 2) 



where represents the increase per unit time of the class (n) due to 
condensation of the vapor on the surface of droplets in class (n - 1) over 
the number, which, owing to evaporation, pass from the class (n) to 
the class (n - 1). 

The rate of change of the number of drops of a given class is deter- 
mined from (A. 1) and (A. 2) to be 



0 f 

T7 = J ” J i 

0 1 n n - 1 



(A. 3) 



Frenkel employs the method of Zeldovich which utilizes the feature of 
the sharp maximum of at n* to solve the equation. By considering 
only values of n > 10 say, quantities appearing in the equation will 
only vary slightly when n is changed by 1 and can be treated as functions 
of a continuous variable. Equation (A. 3) can then be integrated and the 
result is presented as equation (27b) on page 396 of Frenkel (13). 



J = D(n*)| exp 



A4> (n*) 



kT 



I 



2v kT 



(A. 4) 



where - number of molecules in the supersaturated vapor 

E>(n*) = S(n*) j3 - diffusion coefficient = (4 tt )^^ 3^^ V n *^/3 ^ 

p ^ 

/3= £ 

V 27T mkT 

4 2 

A4> = y tt a r* 



1 

7 ~ 3 



(>,. - 0 )/ n* 

liq v ' 



substituting these values into (A. 4) produces the equation 

. 2 ' 

P 



J = 2r* 



kT 



jexp 



4rff r* 

~ 3kT 




i ( V 



0 )■ 
v 



mir 



U 5) 
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or, writing in the form of Frenkel's final equation (28a) 

A 



n*l - or* 2 -E-Jexo - 477 arJ 
N 2 kT\ exp 3kT 
n 



(i ijq 



. - <t> ) 
v 



m 



ir 



(A. 6) 



In comparing (A. 6) with Frenkel's equation (28a), it is noted that Frenkel 
does have an extra factor of n in the coefficient as noted by Courtney (9) 
and others. 



The desired equation for this analysis is obtained by substituting 

2a V.. 

iiq 



the value of - <f>^) from the critical radius equation r* 



and the value n* 



14 3 

— — g 71 r * into equation (A. 5) 



J ■ 2r * 2 (is 



liq 

2 / 



exp - 






2r* 



2r* 



2 V liq 




- 2a V v V.. 3 
1 liq liq 



4tt r’ 



p \“l 4 ttc r * 

kT eXP " 3kT 



<t>,- ~ 

liq v 



J = 



m) V liqVH eXt) - 



47ra r*' 
3kT 



To eliminate the necessity of calculating molecular values of 

*A 



(A. 7) 



V,. = 



, . „ T and m = tF - , the equation for computing J in this 
liq N A P L N ' 



analysis is written 



, _|-P _| 2 J-jlmL 

' kT/ p L U 



l1/2 



exp 



47r<7 (r*)' 
3kT 



(A. 8) 



Appendix B 

IBM 7090 DIGITAL COMPUTER PROGRAM FOR 
CONDENSATION IN HIGHLY EXPANDED FLOWS 

Written in Michigan Algorithmic Decoder (MAD) Language 

This program is valid for a pure vapor with constant specific heat 
and constant latent heat of vaporization. 

The cgs system of units is utilized throughout and cgs values for 
the universal gas constant, Boltzmann’s gas constant, and Avogadro's 
number have been included within the program. 

The symbolic names of parameters used in the program are listed 
and defined at the beginning of each applicable section of the program. 

The units for and input data required for utilization of this program 
are presented in Appendix C. 
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SCOMPIlE MAD .PRINT OBJECT*PUNCH OBJECT MAIN 001MA1N 000 

R AERONAJT IC AL ENGINEERING 

R' NOZZLE CONDENSATION' RESEARCH 

R 

R " MAIN PROGRAM' 

R 

RRESEARCH — j'. GRIFFIN " PROGRAMMI NG — L. HORDING "JULY 1963 

R 

' "PROGRAM COMMON GAMMA ♦ MUT BV'cT'CPV lV~RH0l7'' S IGMA 
R 

START READ" DATA 

R 

RTHlS ^AT A TH^LO~rNCLlJDE 
R 

R 'P ZERO" =''TnTTTal TRES' SURE* 

R RHZERO = INITIAL DENSITY 
R T ZERO '= 1 INI T I A L _ TEMPER AT URE 

R PE = I SENTROP I C EXIT PRESSURE. THIS IS USED AS AN 

~ R INITIAL APPROXIMATION TO THE SATURATION PRESSURE. 

_R_ 

R GAMMA" '«* RATIO OF' THE SPECIF I C 'HEA f S 
R MU = MOLECULAR WEIGHT OF THE VAPOR 

R" CP = SPECTFTC HEAT AT '“C O'NST A Nf"~ PRESSURE 

R L = LATENT HEAT (NOTE VAPOR PARAMETERS SUBROUTINE) 

FT SIGMA = SURF A"C E“T ETNSTOTTTNOTE VAPORS PARAMETERS SUBROUTINE ) 

R 

R" C * **SAT UR A T I ON C UR V E " CONS T a NT 

R B = SATURATION CURVE CONSTANT 

R ALPHA’” = ACCOMODATION COEFFICIENT ~~ ~ ' 

R 

'R DELX =‘ INCREMENTAL'“STEP T'N^X. THE NOZZLE 7 OWfn ONS "' W T"LL 
R BE COMPUTED AT XCON+DELX , XC0N+2^DELX_» . . . WHERE 

R ~ XCON IS* "THE CONDENSATION POINT”. 

R X RANGE = LENGTH OF INTERVAL. STARTING AT THE CONDENSATION 

“R POINT." IN WHICH THE FLOW CONDITIONS WILL~'BE~“' 

R COMPUTED. 

R “XFOrNT = THE’ CONDENSATION POINT" IS ASSIGNED A SUBSCR I PT "OF 
R ZERO. ONLY THOSE POINTS WITH SUBSCRIPTS XPOINT. 

K 2*XP0 1'NT > . . . ’WI LL BE PRINTED." ' 

R 

"R T RANGE" = P R"0G R AM""W I L L LOOK FOR A CONDENSATION POINT IN THE" 

R INTERVAL ( TSAT .TSAT-TRANGE ) . WHERE TSAT IS THE 

R COMPETED SATURATION TEMPERATURE. 

_R DJLTAT = INITIAL __S T E P SIZE USED IN LOOK I NG_FOR CONDENSATION 

r "point in the Tnte'rval '(Tsat"»Tsat-trang"e 

R PERCNT = IF EPSLON(T) EXCEEDS PERCNT IT IS ASSUMED THAT 
R ' CONDENSATION HAS STARTED. ' ' 

R 

PRINT COMMENT $1 " ' '' NOZ 

_1ZLE_ CONDENSATION RESEARCH $ 

"" P'RIN't'"COMMENT' $6" ' "INITIAL CONDI TIONSS 

PRINT RESULTS P ZE RO ♦ RHZ ERO . TZ ERO , P E 

""PRINT' COMMENTED" “ CONSTANTS' DEPENDING ON THE VAPORS 

PRINT RESULTS MU »CP . GAMMA , S I GMA . L 

PRINT COMMENT "TO SATURATION CURVE 'CONST ANTS $ 

PRINT RESULTS C.B 

prinT'TOmmEnt ''$0 PTTogram''paRameterss‘ “ 

PRINT RESULTS DELX , XRANGE » XPOI N T . ALPHA 
ETXTCuTE'NOZ ZEE . ( $THR OATS' ."OTTA'S T A R ) " 
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R FIND SATURATION PUINI BY MATCHING ISENTROPIC PRESSURE 

R SOLUTION Tu T M t SATURATION CURVE PRESSURE. 

k 

P = P E 

Z= GAMMA/ ( GAMMA- 1 . ) 

El T=C/ (jB-ELOG. (P )) 

PSAT =PZFRO* ( (T/TZERO) . P. Z ) 

TEST=.ABS.( 1.-PSAT/PJ-.001 
WHENFVFR TEST .L.O. » I RANmFER 10 ^2 
P=P+.5*(PSAT-P) 

TRANSFER TO SI 
S 2 TSAT = T 

EXECUTE IFLOWS. (ASTAR.CP .GAMMA .PZERO. RHZERO.TZERO) 

EXECUTE I FLOW. (ASAT.MSAT .PS AT .RHOSAT . TSAT. US AT ) 

Z = MSAT 

EXECUTE NOZZLE .($ I NVERSS ,Z .ASA T ) 

WHENEVER Z ( 1 ). E . 1 .. TRANSFER TO S3 

PRINT COMMENT $4 THE SATURATION POINT FOUND IS NOT INSIDE 

1 THE NOZZLE. THE PROGRAM WILL CONTINUE. S 
S3' X S A T = Z 

MDOT=ASAT*RHOSAT*USAT 

PRINT COMMENT 3> 2 COMPUTED SATURATION CONDITIONSS 

PRINT RESULTS PSA T » RHOSAT . TSAT 

PRINT RESULTS MSAT.USAT 
PRINT RESULTS XSAT.ASAT 
PRINT RESULTS MDCT 

R 

R FIND CONDENSATION POINT BY MATCH ING ISENTROPIC FLOW 

R SOLUTION TO THE NUCLEAT ION THEORY SOLUTION GIVING THE 

R AMOUNT OF CONDENSATE AS A FUNCTION OF AN ASSUMED FLOW. 

DELT =DELTAT 

I TERF=0. 

TMIN=TSAT-TRANGE 

WHENFVFR TM I N . L . 0 . » TM. I N= 0 . 

P"RINT COMMENT 42 PARAMETERS DETERMINING CONDENSATIO 

IN POINT SEARCHi 

FETNT RESULTS TSAT .TMIN.DLLT 

PRINT RESULTS PERC.NT . 

“R 

R 

34 T=TSAf' 

I T ER =0 . „ 

T = T-DELT 

WHENEVER T .L. TM I N 

“ Z= 1 . — 

S 6 ITERF=ITERF + 1. 

— DTlT=.5*0FLT 

WHENEVER ITERF . L . 1 0 « » TR A N oFER Tu S 4 

WHENEVER Z . E . 0 .♦ TRANSFER TO S7 

PRINT COMMENT 40 TEN PASSES OVER T EMPE RATURE RANGE SPECIFIED. 

1 NO CONDENSATION POINT INDICATED. S 

TRANSFFR TO START 

END OF CONDITIONAL 
EXECUTF I FLOW. (A .M.P.RHO.T »U) 

Z = M 

EXECUTE NOZZLE. < SINVEhLii.Z.A) 

WHENEVER Z ( 1 ) . E . 0. 



-- 1 4 -_ 



z=o. 

TRANSFER TO S6 

s7 prin't comment so 'ten "passes "over The full range' of the’ nozzle 

1_* NO CONDENSATION POINT INDICATED. S 

'TRANSFER TO START' "" 

END OF CONDI TIONAL 

EXECUTE NOZZLE. ( $ ARE A $ , X-DELX ♦ Z ) 

DElA = A-Z ' 

R 

L = FL • ( T ) 

TDROP=C/ (B-ELOG. (P) ) 

RHOL=FRH0L. ( TDROP ) ' — 

SIGMA = FSIGMA. ( TDROP) 

EXECUTE NUCLEI.' (A, DELX ."MOOT » P ♦ T » RAD » NDOT ♦ DELG)' 

R 

R ' COMPUTATION OF 7 ' EPS I LON ( T) 

R 

' E P S L 0 N ='( GAMMA- ( 1 ./ (M*M) ) )7TgAMMA- l'.1 

WH ENEV ER M. L. 1 . ♦ E PSLON=l . 

EPSLON = A»f 1 .'- ( L/ ("CP*f ) 7*EP"SLON ) 

EPSLON=EPSLON* (DELG /DEL A) 

EPSl6n= . ABS • EPSLON 
R 

whenever eP'slon. l.peRcnt .transfer to ~s~5 ~ 

t=t+delt 

DELT = .5*DELT 

ITER-ITER+1. 

WHENEVER I TER. L. LI TER .TRANSFER TO S5 
R 

R LIT ETThALF-INT E R V AT" STEPS'" HA V E"B"E E N~ PER F oTTmED - AND 

R A CONDENSATION TEMPERATURE HAS BEEN FOUND. 

R 

PRINT COMMENT $0 HAVE FOUND TEMPERATURE THAT SATISFIES THE CO 
1NDIT IONS' FOR THE CONDENSATION POINTS 
PRINT RESULTS T.DELT 
TCON=T 

EXECUTE IFLOW. < ACON ,MCON .PCON , RHOCON . TCON . UCON ) 

Z = MCON - " 

EXECUTE NOZZLE. (SINVERSS.Z. ACON) 

WHEREVER - ZTl j.E.l., TRANSFER TO 'S8 

PRINT COMMENT SA THE CONDENSATION POINT FOUND IS NOT INSI 

IDE THE NOZZLE. ‘THE PROGRAM WILL CONTINUE. $ 

58 XCON =Z 

L=FL.(T) 

TDROP=C/ (B-ELOG. (PCON) ) 

RHOL=ERHOL“; ( TDROP ) 

S I GMA = FS I GM A • ( TDROP )_ 

~ execute - Nuclei'. ( acon .delx.mdoT .pcon .tcon.rad .ndot.delg j 

G= DELG 

FRT'nT COMMENT '$ A ~ COMPUTED 'CONDENSATION CONDITIONS* 

PRINT RESULTS PCON.RHOCON. TCON 
- PEINT" RESULT'S MCON .UCON 
PRINT RESULTS XCON. ACON 
PRINT RESULTS S I GMA . RhOL . L , CP 

PRINT COMMENT $_NUCLEATION THEORY QUANT I T I ESS 
‘TR InT "RESULTS' RAD .NDOT '.DELG. G 
R 

“RlNlTIALIZE - FOR" 'NOZZLE COMPUTATIONS STARTING AT THE CONDENS- 
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>" 2 ¥ 



RATION POINT . THE VECTOR X CONTAINS THE POINTS ALONG THE 
RAXIS OF THE NOZZLE AT WHICH THE QUANTITIES DESIRED ARE TO 
RbE COMPUTED. NOTE THAT IT Is NOT NECESSARY THAT THE SPACING 
UNIFORM SINCE A VECTOR DELX IS PROVIDED. AT EACH POINT 
THE MESH WE WANT TO CALCULATE 



RbE 

ROF 

R 



P 

RHQ 

T 

U 



DELG 

G 



R 
R 
R 
R 
R 
R 
R 

ROlLR ( I ) 
R 

R RAD ( I ) 
RNDOT ( I ) 
R 



THE LA^T MESH 
VAPOR THAT IS 



THE 



LAST 
I 



POINT 
IN THE 

POINT * OF 



X ( I ) 

ARE NEWLY 



VAPOR PRESSURE 
MIXTURE DENSITY 

temperature 
flow velocity 

MACH NUMBER 

CONDENSATE INCREASE oINCE 
FRACTION* BY MASS* OF THE 
CONDENSED sTATE 
THE INCREASE IN RADIUS* SINCE 
THE DROPS INITIALLY FORMED AT POINT 
RADIUS OF DROPS INITIALLY FORMED AT 
THE NUMBER OF DROPS PER SECOND THAT 
FORMED FROM X( 1-1 ) TO X ( I ) 

SURFACE TENSION 
LIQU ID DENSITY 
SPECIFIC HEAT 
LATENT HEAT 

THE IDENTIFYING QUANTITIES 
X = POINT OF EVALUATION 
D"! STANCE FROM 
NOZZLE AREA 

DELA = AREA INCREASE FROM LAST MESH POINT 

THE POINTS* THEIR INCREMENTS AND THE INTEGER INDEX N 
KEPT IN CORE* ALL THE REST OF "THE I N~F 0 R M AT i 0 N I S " MAIN- 
RTAINED ON TAPE. ALL QUANTITIES FOR POINT N* HOWEVER* ARE 

'obtain Those for point n'+i* hence They' are 

IN CORE AS UNINDEXED QUANTITIES WITH THE ABOVE 



SIGMA 

RHOL 

CP 

L 



ralso 

R 

R “DELX“^ 
R A = 



AT CONSTANT PRESSURE 



LAST MESH'POINT 



R 

RONLY 

RARE 



RREQU I RED TO 
RTEMPORAR I LY 
RNAMEsV 



R INITIALIZE TO CONDENSATION POINT 

R 

? T = TCON 
P = PCON 

rho=rhocon“ 



U = UCON 
M=MCON 
X = XC ON 
A = AC ON 

EXECUTE NOZZLE. ( SARE AS * XCON-D ELX * A ( 1) ) 
D*E LA = A- A (1 1 
N = 0 

PPC I N T = 6 



EX ECU TF IOCTRL. 

“DElu=.oo6i*ucon 

XL IM I_T = X C 0 N + XR A N G E 

R“ 

R 

U( 1 ) =u+deUu — 

WHEN E VFR U( 1 ) .L.O. 

“Print comment' s4 Th^TnitTal “approximation to the velo 

1CITY AT THE NEXT MESH POINT IS NEGATIVE. S 
TRANSFER TO S 34 
END OF CONDITIONAL 

■R 
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S 30 



331 



S3T 



S33 



R BEGIN iteration on delg until nucleatign theory and 

R FLOW EQUATIONS AGREE. A MAXIMUM OF GI1LR ITERATIONS ARE 

R ALLOWED PER MESH POINT. 

R 



I TFR=GI TER 



N = N+ 1 

WHENFVFR N.G.1C00 

PRINT COMMENT $4 NUMBER OF MESH POINTS LXC.FE DS THE PRO 

1GRAVMED LIMIT.5 

TRANSFER TO S 3 4 
END OF COND I T I ON AL 
D E L X ( N ) = D F L X 
X ( N ) = X t M-l ) +DELX ( N ) 

WHENEVER X.G.XLlMlT 

PRINT COMMENT 2>4 THE NEXT POINT WOULD EXCEED THE RwNUr. 

1 DESIRED. 5 

Transfer to S34 

END OF COND I T ICNAL 
EXECUTE NOZZLE. ( $AREAi>*X(N) , A ( 1 ) ) 

WHENEVER At 1 ) .F.O. 

PRINT COMMENT $4 NOZZLE SUBROUTINE ERROR INDICATION - 

1KESH POINT DOES NOT LIE INSIDE NOZZLE. $ 

PRINT RESULTS X(N> 

TRANSFER TO S34 
END OF CONDITIONAL 
DELAt 1 ) = A ( 1 ) - A 
DELG ( 1 ) =DELG 
G ( 1 ) =G+DELG ( 1 ) 

EXECUTE ' CFLOW. < A ( 1 > .DEL A ( 1 ) ,G < 1 ) . DELG ( 1 ) » P » RHO , T » U *M ) 



L = FL . ( T ( 1 ) ) 

TDR0P = C/ (B-ELOG. ( P ( .1 ) ) ) 
RHOL=FRHOL . ( TDROP ) 



S I GMA = FS I GMA . ( TDROP ) 

TEST =G ( 1 ) 

Execute nucle 2. < x t n > ,delx < ni , a m . delai i i .p t i i »rho< i i »t t i i 

1U( 1 ) »M( 1 ) .MOOT .ALPHA »N . RAD , DR AD » NDOT . DLLG ( 1 ) ) 

“ C( 1 ) =G+DELG( 1 ) 



TEST =. ASS. ( TEST-Gt 1 ) ) 

WHENEVER TEST . L . .00001 . TRANSFER TO S32 

I T ER = I T F R- 1 . 



WHENEVER I TER. GF.O. .TRANSFER TO S31 

PRINT COMMENT $ S 

_ X=X ( N) 

DELX =DFL X ( N ) 

A = A ( 1 ) 

DELA = DFLA( 1 ) 

P = P( 1 ) 

RHO = RHC ( i ) 



T = T ( 1 ) 



DE LU = U ( 1 ) - U 
U = U ( 1 ) 



M = M ( 1 ) 

DELG = DFLG< 1 ) 

G=G ( 1 ) 

“THROUGH S^S.FOR I = 0 . 1 . I . G . ( N- 1 ) 
RAD ( I ) =RAD ( I ) +DRAD ( I ) 

EXECUTE IOCTRL. 

"FT 



or or a: !(v o: 



534 
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transfer TO S29 
R 

transfer to start 



INTERNAL FUNCTION FOR OUTPUT OF THE FLOW QUANTITIES AT 
MFSH POINT X. LATER THIS SECTION MAY ALSO SAVETHIS 
INFORMATION FOR FURTHER PROCESSING. 

R ' 

INTERNAL FUNCTION 
R 

ENJRY TO IOCTRL. 

R 

WHENEVER N.E.PPOINT 

PRINT COMMENT $4 ‘ ' NOZZLE CONDITIONS — S 

PRINT RESULTS N ♦ X * DELX , A . DELA , P , RhO . T . U *M * RAD < N ) ,NDOT < N ) » 

IDE LG »G 

WHENEVER N.G.O 

PR IN T RESULTS RADIO )...RAD(N-1) 

_ END OF CONDITIONAL 

PPOI NT =N + XPO I NT 

END_OF CONDITIONAL 

FUNCTION RETURN 
R 

~END'OF FUNCTION 
R_ 

INTEGER I ,N .PPOI NT 

DIMENSION DELX ( 1000 ) > DR AD ( 1 000 ) »NDOT ( 1 000 ) j.R AD ( 1 000 ) » X ( 1 0 00 ) 
DIMENSION' All') , DELA I 1 ) * DELG I 1 ) ,G( 1 ) .MI 1 ) * PI 1 i'.RHOI 1 ) * T I 1 ) , 

1 U I 1 ) » Z I 1 ) 

VECTOR VALUES XPOINT = 1 
VECTOR VALUES LITER = 10. 

' VECTOR ' VALUES G I TER = 10. 

R 

END’ OF PROGRAM 
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IFLOWOOl IFLOWOO 



{.COMPILE MAD * PRINT OBJECT .PUNCH OBJECT 

R __ AERONAUTICAL ENGINEERING 

' R~ " "" "NOZZLE CONDENSATION RESEARCH 
R 

R ‘ I SENTROP I C FLOW 

R 

RRESEARCH — J. GRIFFIN PROGRAMMING-- L. "HARDING JULY 1963 

R 

RlSENTROPIC FLOW EQUATIONS FOR COMPUT f N G 'P R Fs SUR i7 " M A C H NUMBER 
RAREA, VELOCITY AND DENSITY GIVEN THE TEMPERATURE AND THE 

RFOLLOWING PARAMETERS 

R ASTAR - THROAT AREA OF NOZZLE 

~R - CP ^ SPECIFIC HEAT AT CONSTANT PRESSURE ~~~ 

R GAMMA - VAPOR EXPONENT 

R’ PZERO“- INITIAL PRESSURE ‘ 

R RHZERO - INITIAL DENSITY 

"R tzero - initial temperature" 

R 

EXTERNAL FUNCT ION" ( AlTAT7A3'»A4TA5»A6 ) " 

R 

"R “THIS ENTRY PROVIDES 'FOR PICKING UP THE PARAMETERS LISTED 
R ABOVE. THEY ARE THEN SAVED INTERNALLY FOR LATER USE. 

R" ' ‘ 

ENTRY TO IFLOWS. 

■~A'$tar*ai — 

CP = A2 
GA’MMA = A3 
PZER0=A4 

RHZERO = A5 

TZERO=A6 
FUNCTION RETURN 
R 

TT 

R 
R 
R 
R 
R 
R 
R 
"R 
R 
R 

ENTRY TO IFLOW. 

T = A5 ’ ~ ...... .... 

R PRESSURE 

— EXP = GAMMA/ ( GAMMA-1 . ) 

Z=(T/TZERO) .P. EXP 
~'A'3 = PZERO*Z 

R VELOCITY AND MACH NUMBER 

“ Z'=2T*TTZERO’-T ) " ' 

A6=SQRT • (CP*Z ) 

~T2Z7TGAMMA-T.it ' 

A2=SQRT . ( Z/T ) 

~R D E N S I T Y - ATnT'NOZ Z L'E 'ARE TT - 

EXP= 1 ./ ( GAMMA-1 . ) 

Z = 1 . + . 5 : ® r A'2 r F'A77TXR 
A4=RHZER0/(Z .P. EXP) 

— ETF=“r5^T GAMMA+TTT*EXR 



THIS ENTRY COMPUTES THE UNKNOWN QUANTITIES 
A 1 = NOZZLE AREA 

A 2 = MACH NUMBER 
A3 = PRESSURE 
A 4 = ' DENSTTYT ~RHO T~ 

A6 = VELOCITY 

IN TERMS OF THE PARAMETERS SAVED "FROM THE LAST CALL 

OF 'IFLOWS' AND 

A5 = TEMPERATURE 

NOTE THAT THE LIST OF ARGUMENTS IS ALPHABETICAL. 
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Z = 2 . * 17 ( GAMMA+ 1 . ) 
_Z=(Z .P. EXP ) /A2 

AT=astar*z 

R 

“FUNCTION (RETURN 
R 

END OF FUNCTION 
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iCOMP I L E MAD, PRINT OBJECT, PUNCH OBJECT CFLOWOOICFLOwOO* 

R AERONAUTICAL ENGINEtkING 

R NOZZLE CONDENSATION RESEARCH 

R 

R CONDENSATION FLOW 

R 

RRESEARCH-- J, GRIFFIN PROGRAMMING — L. HARDING JULY 1963 

R 

RCOMPU T A T I ON OF VALUES OF PRESSURE ,DENSI TY , TEMPERATURE , MACH 
RNUMOER AND VFLCCITY THAT SATISFY THE FLOW EQUATIONS THAT 
RAPPLY AFTER THE ONSET OF CONDENSATION. IT IS ASSUMED THAT WF 
RHAVE COMPUTLD THE CONDITIONS AT X AND DESIRE THOSE AT X+DELX 

R the arguments required are 



R 


DELA 


- 


AREA INCREASE FROM X TO X+DELX 




R 


A 


= 


AREA AT X+DELX 




R 


dflg 


= 


CONDENSATE MASS INCREASE FROM X TO 


X+DELX 


R 


G 


- 


CONDENSATE MASS AT X+DELX 




R 


P 


- 


PRESSURE AT X 




R 


RHC 


= 


DENSITY AT X 




R 


T 


= 


TEMPERATURE AT X 




R 


U 


= 


VELOCITY AT X 




R 


M 


= 


MACH NUMBER AT X 




RTHE 


PARAMETERS CP*L AND MU ARE REQUIRED FROM 


PROGRAM COMMON 



RTHE FUNCTION FL.(T) MUST BE USED HOWEVER SINCE ITERATION ON 
RT IS NECESSARY. THE UNKNOWN QUANTITIES ARE RETURNED IN P(l), 
RRHC ( 1 ) ,T ( 1 ) »U< 1 ) AND M(l). IT IS ASSUMED THAT U ( 1 ) CONTAINS 
RAN INITIAL APPROXIMATION TO THE VELOCITY AT THt NEXT MESH 
RPOINT ON ENTRY. 

R 

EXTERNAL FUNCTION ( A , DELA , G , DELG , P A , RHOA , T A , UA * MA ) 

R 

PROGRAM COMMON GAMMA, MU, B, C, CP, L, RHOL, SIGMA 
R 

ENTRY TO CFLOW. 

U=UA ( 1 ) 

DELU=U-UA 
MDGG=DELG/ ( 1 . -G ) 

DQA=DELA/A 
L = FL . ( TA ) 

R_ 

“THROUGH S2 > FOR I = 1 , 1 , I . G . 1 00 
DQU=DELU/U 
R 

R COMPUTATION OF DCLT AND T FROM THE E N t R G Y EQUATION 

'R ' 

DELT=(L*DELG-U*DFLU) /CP 
T = T A+DFLT 
L=FL . ( T ) 

“R 

R COMPUTE DELRHO/RHO FROM CONTINUITY EQUATION 

DQRHO=DQA+DQU 

UQRHO=-DGRHO 

R 

~ R COMPUTE DELT7P FRUM“ EOUAT I ON OF SI ATE 

R _ 

DTTP^ DTJRh'O+TdEl T/ T 

R 

— R WE CAN NOW' COMPUTE THE VELOCITY CORRESPONDING TO 



R THESE DIFFERENCE QUOTIENTS FROM THE MOMENTUM EQUATION 

R 

DELU=-(DQP*R*T*( 1.-6) ) / (MU*U) 

NEWU=UA+DELU 
TEST=.ABS. ( l.-NEWU/U) 

U=NEWU 

WHENEVER TEST • L. EPSLON * TRANSFER TO S3 

R 

R ITERATION LIMIT HAS BEEN REACHED WITHOUT CONVERGENCE 

R 

PRINT COMMENT $ 'CFLOW' - ITERATION FAILED TO CONVERGES 
EXECUTE ERROR. 

T? 

R CONVERGENCE CRITERIA SATISFIED 

R 

TA ( 1 ) = T 
UA ( 1 )=U 

DQRHQ=1.-DQRH0 

DQP= 1 . -DQP 

PA ( 1 ) =PA/DQP 

RHOA ( 1 ) =RHOA/DQRHO 

NEWT =R»GAMMA*T/MU 

MA< 1 )=U/SQRT. (NEWT > 

FUNCTION RETURN 
R 

VECTOR VALUES R = 8 . 3 1 4E + 07 

VECTOR VALUES E PSLON = 1 . E-06 
INTEGER I 

R 

END OF FUNCTION 
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iCOMPlLE MAD, PRINT OBJECT, PUNCH OBJECT NUCL 1 00 1 NUCL 1 00* 

R AERONAUTICAL ENGINEERING 

R NOZZLE CONDENSATION RESEARCH 

R 

R NUCL E A T I ON THEORY - i 

R 

RRESEARCH — J. GRIFFIN PROGRAMMING — L. HARDING JULY 1 9o3 

R 

RCOMPu T AT ION OF THE MASS DUE TO INITIAL CONDENSATION AS GIVEN 
RtSY NUCLE AT I ON THEORY. THE HOLLOWING QUANTITIES MUST bE AVAIL- 
RABLE. 

R ( 1 ) CONSTANTS 

R NA = AVOGADROS NUMBER (MOLECULES/GRAM MOLE) 

R K = BOLTZMANN CONSTANT (DYNE*CM/DEGREEK) 

R R = UNIVERSAL GAS CONSTANT ( DYNE*CM/ GRAM MOLE*DEGREEK. ) 

R ( 2 ) PARAMETERS 

R B SATURATION CURVE CONSTANT 

R C SATURATION CURVE CONSTANT 

R DEL X = DELTA X FOR VOLUME 

R MOOT = RHOSAT*ASAT*USAT 

R MU = MOLECULAR WEIGHT OF VAPOR 

R<3) VARIABLES 
K A = NOZZLE AREA 

R P = PRESSURE 

R T = TEMPERATURE 

R RHCL = LIQUID DENSITY DEPENDENT ON TEMPERATURE 
R SIGMA = SURFACE TENSION DEPENDENT ON TEMPERATURE 
R 

EXTERNAL FUNCTION ( A 1 ,A2 . A3 . A4 , AS , A6 , A7 ♦ AS ) 

R 

PROGRAM COMMON GAMMA. MU. 6. C. CP. L. RHOL . gIGMA 
R 

RCOMPUTE THE RADIUS OF THE DROPS CONDENSING UNDER THESC 
RCCNDITIONS. THEIR NUMBER AND HENCE THE CONDENSATE MASS 
RFORMBD IN A VOLUME ELEMENT OF LENGTH DELX. THE PROGRAM 
RCOMMON PARAMETERS RHCL AND SIGMA ARE ASSUMED TO CORRESPOND 
RTO THE DROP TEMPERATURE FOR THE GIVEN FLOW CONDITIONS. 

R 

ENTRY TO NUCLEI. 

A = A 1 
DELX = A2 
MOOT =A3 
P = A4 
T = A5 
R " 

R COMPUTE RADIUS < CM. ) OF THE DROPS THAT ARE CONDENSING 

R 

RADIUS 1 ( 2.*SIGMA*MU ) / ( RHOL* T ) 

RADIUS = RAD I US/ ( R* ( ELOG. ( P)-B + C/T ) ) 

R 

ft" COMPUTE THE NUMBER OF DROPS OF THIS SIZE THAT ARE 
R CONDENSING PER CENTIMETER CUbED PER SECOND. 

R 

TtXP=RADIUS*RADIUS/K 

TEXP=- (4. 1887902*5 I GMA*TEXP )/f' 

N=SQRT. ( SIGMA*VU/NA( 1) ) 

N=ft*N7TK ( 1 ) *T f ’ ‘ ' ' ~ " 

N=P*N*EXP. ( TEXP ) / ( K*T*RHOL ) 

N = • 7 9 78 8 46 *N 
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ndot=n*a*delx 

R 

R COMPUTATION 'OF PERCENT OF LIQUID MASS. DELTAG 

R 

DELG=4. 1 887902*RHOL*NDOT*RADIUS 
DELG=bELG*RADI US*RADIUS/MDOT 
R 

_A_6f RADIUS 

A 7 = N DO T 

A8= DELG 

FUNCTION RETURN 
R 

VECTOR VALUES NA =6 . 02 7E+2 3 . 6 . 02 7E+03 

VECTOR VALUES, K= 1 . 379E- 1 6 . 1 * 37 9 E-06 

vector Values ‘r = 8.314E+07 

R 

END OF FUNCTION 
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SCOMPILE MAD»PR’lNT OBJECT.PUNCH OBJECT NUCL200 1NUCL200* 

_R AERONAUT ICA L ENGINEERING 



R 

R 


NOZZLE CONDENSATION RESEARCH 




R 

R 


NUCLEAT I ON THEORY - 2 




rResearch — 


J. GRIFFIN PROGRAMMING — L. HARDING 


JULY 1963 



R 

RCOMPUTAT ION O F" T H E Q U A N T f T I E S ARISING DUE TO CONDENSATION 



RUNDER THE 


FOLLOWING CONDITIONS 




R 


X 


= 


NOZZLE COORDINATE 




R 


DELX 


= 


INCREMENT SINCE LAST MESH POINT 




R 


A 


= 


NOZZLE AREA AT X 




R 


DELA 


= 


NOZZLE AREA INCREASE SINCE X-DELX 




R 


P 


= 


PRESSURE AT X 




R 


RHO 


= 


MIXTURE DENSITY AT X 




R 


T 


= 


TEMPERATURE AT X 




R 


U 


= 


VELOCITY OF FLOW AT X 




R 


M 


= 


MACH NUMBER AT X 




R 


MDOT 


= 


TOTAL MASS FLOW RATE 




R 


ALPHA 


= 


ACCOMODATION COEFFICIENT 




R 


N 


= 


MESH POINT NUMBER SINCE CONDENSATION. THE POINT 





R OF CONDENSATION IS MESH POINT ZERO. 

R THE QUANTITIES TO BE COMPUTED ARE 



"R RAD(N) = RADIUS OF NEW DROPS FORMED AT X. ( RADI I } IS RADIUS 
R AT X-DELX OF DROPS INITIALLY FORMED AT MESH POINT 

R ' I AND IS NECESSARY INPUT ) 

RNDOT(N) = NUMBER OF DROPS OF RADIUS RADI N ) FORMED BETWEEN 
R X-DELX AND X. ( NDOT ( I ) IS THE NUMBER OF DROPS 

R INITIALLY FORMED BETWEEN MESH POINTS 1-1 AND I AND 

R IS NECESSARY INPUT. ) 

RDELR ( 0) _= DELRI I ) IS THE RADI US IN CR E AS E BETWEEN X-DELX AND X 
R THRU “ "OF fHOSE DROPS 'INITIALLY FORMED BETWEEN MESH 
RDELR(N-l) POINTS 1-1 AND I» THESE QUANTITIES MUST ALL BE 
R COMPUTED. ) 

R DELG = CONDENSATE MA S S INCREASE FROM X-DELX TO X DUE TO 
R NEWLY FORMED DROPS AND GROWTH OF OLD DROPS. 

R 

RTHE ARGUMENTS ARE THE QUANTITIES GIVEN ABOVE IN THAT ORDER. 
RNOTE THAT THE ZEROTH WORD OF THE VECTOR ARGUMENTS ARE TO BE 
RGIVEN. FURTHER THESE VECTORS CONTAIN REQUIRED INPUT AS WELL 
RAS SPACE TO PUT THE DESIRED OUTPUT. FURTHER THE USUAL 

rprogram common parameters are assumed to correspond to the 

RARGUMENT TEMPERATURE. 

R 

EXTERNAL FUNCTION ( X . DELX » A ♦ DELA , P » RHO » T , U »M .MDOT . ALPHA , N . 

1 RAD .DELR .NDOT .DELG ) 

R 

PROGRAM COMMON GAMMA. MU. B. C. CP. L. RHOL. SIGMA 
INTEGER N 

"R 

ENTRY TO NUCLE2 . 

"R 

EXECUTE NUCLEI. ( A .DELX .MDOT ♦ P . T . RAD ( N ). NDOT ( N ). DELG ) 

R " 

R WE HAVE NOW ONLY TO COMPUTE THE RADIAL INCREMENTS AND 

"K THEIR CONTRIBUTION TO THE CONDENSATE MASS INCREASE. THE 

R DROP TEMPERATURE IS TAKEN AS THE SATURATION CURVE 

~R TEMPERATURE CORRESPONDING TO THE GIVEN PRESSURE. 
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R 

_Tl) R CP = C/ (D-ELOG. ( P) ) 

DE L R = 79 78 • 8 46 * SORT «" ( 1 . 379* . 60277 (T*MU )) * ( foROP-T 7 
DE LR =DEL X* ( ALPHA/ L )* ( P/RHOL ) * ( DELR/U ) 

THROUGH SI, FOR I = 1 * 1 » I .G. ( N-l ) 

OELR ( 1 ) =DELR 



COMPUTATION OF RESULTANT TOTAL CONDENSATE MASS INCREASE 

R ~ ' 

DELG2=0. 

THROUGH S2'» FOR I =0 , 1 » I . G. ( N-l ) 

DELG2 = DELG2+ R AD ( I )*RAD( I ) *NDOT ( I ) *DELR ( I ) 

' DELG=DELG+12.5663706*RHOL*DELG2/MDOT 

R 

’ FUNCTION RETURN 
R 

INTEGER I 
R 



END OF FUNCTION 
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BCuMPILL 



>'AD, PRINT OBJECT , P U N C H OBJECT NOZLEOOlNGZLtOO* 

R AERONAUTICAL ENGINEERING 

R NOZZLE CONDENSATION RESEARCH 

R 

R NOZZLE 

R 



R RE SEARCH — J. GRIFFIN PROGRAMMING — L. HARDING JULY 1963 

R 

RORIGIN OF THE COORDINATE SYSTEM FOR THE NOZZLE IS THE CENTFR 
ROF THE THROAT * I.E.f A ( 0 ) = ASTAR • THE VARIABLE X INCREASES 

RPOSITIVELY IN THE DIRECTION OF FLOW AND IS NEGATIVF UPSTREAM 
RFROM THE THROAT. THE NOZZLE PARAMETERS ARE 



R 

R 

R 

R 

R 

RTHF 

R 



ASTAR = THROAT AREA 
XMIN = X-COORDINATE OF INTAKE 
XM AX = X-CO ORDINATE OF EXIT 
INANG = INTAKE HALF ANGLE IN DEGREES 
OUTAMG = EXIT SIDE HALF ANGLE IN DEGREES 
ENTRY THROAT RE ADS AND PRINTS THESF PARAMETERS. 



EXTERNAL FUNCTION ( A1 ,A2 ,A3 ) 

R 

ENTRY TO NOZZLE. 

R 

WHFN EVER A1 .E. STHROATS 
R 

R RETURN THROAT AREA AND INITIALIZE IF NECCESSARY 

READ DATA 

PRINT COMMENT $0 NOZZLE PAR AME TERSS 

PRINT RESULTS ASTAR* INANG,OUTANG 
PRINT RESULTS XMINfXMAX 
A= I NANG/RADI AN 
I N T A M = S I N • ( A ) /COS. ( A ) 

A = 0 ! J TANG/RADIAN 
OUTTAN=S IN. ( A ) /COS. ( A ) 

A3=ASTAR 
FUNCTION RETURN 
R 

OR WHENEVER A1 .E. SAREAi 
R 

R COMPUTE THE AREA AT A2 AND RETURN IN A3 

WHENEVER (XMIN.L.A2) .AND. (A2.L.0.) 

A= I N T AN* ( . APS. A2 ) 

OR WHENEVER (A2.G.0.) .AND. (A2.L.XMAX) 

A=OUTTAN*A2 
OTHERWISE 
A 3 = 0 • 

WHENEVER A2.E.0*,A3=ASTAR 
FUNCTION RETURN 
END OF CONDITIONAL 
A3=ASTAR+2 • *A 
FUNCTION RETURN 
R 

OR WHENEVER A1 .E. SINVERSS 
R 

R COMPUTE THAT POINT ON THE X-AXIS OF THE NOZZLE THAT 

R HAS AREA * A 3 1 • THIS POINT DEPENDS ON THE MACH NUMBER 

R WHICH MUST BE IN A2 ON ENTRY, SEE THE PROGRAM FOR 

R TnE PARTICULARS. IF A 2 ( 1 ) IS 1. ON RETURN THEN THE 

R SOLUTION IS IN A2 , IF IT IS ZERO THEN THE SOLUTION IS 



R STORED IN A2 BUT DOES NOT LIE 

A= .5 »( A3-A5TAR ) 

WHENEVER A • L • 0 • » A = 0 « ' 

X( 1 ) =1. 

MACH=A2 

WHENEVER MACH.L.l. 

X=-A/ INTAN 

WHENEVER XM I N . G . X » X ( 1 ) = 0 . 
“'OTHERWISE' 

x=a/outtan 

WHENEVER X . G . XMAX » X (1) = 0 . 

END OF CONDITIONAL 
A2 = X 

A2( 1 ) = X < 1 ) 

FUNCTION' RETURN 

R 



R 



END OF CONDITIONAL 



INTEGER A 1 
DIMENSION X ( 1 ) 

VECTOR * V A L U E SR AD I A N = 5 7 . 2 9 5 7 7 9 5 
R 



END OF FUNCTION 



INSIDE THE NOZZLE. 
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SCOMPILE MAD, PRINT OBJECT, PUNCH OBJECT VAPOROO 1VAPOROO* 

R AERONAUTICAL ENGINEERING 

R NOZZLE CONDENSATION RESEARCH ~ 

R 

R '“VAPOR PARAMETERS " 

R 

RRESEARCH — J. GRIFFIN PROGRAMMING — L. HARDING JULY 1963 

R 

RCOMPUTE AS A FUNCTION OF TEMPERATURE I’hTTUnCTToNS 
R RHOL = LIQUID DENSITY 

R L = LATENT HEAT 

R SIGMA = SURFACE TENSION 

RNOTE THAT EACH OF THE QUANTITIES MAY BE COMPUTED SINGLY AND 

R THAT THEIR VALUES ARE STORED IN LOW CORE. TH IS PE R MITS THEM 

rt'o be held constant wTthout''‘in'i t i aTTzTng 'these SUBROUTINES. 

RANY ONE, OR ALL, OF THESE PARAMETERS MAY BE HELD CONSTANT BY 

"RSIMPLY returning the value 'from program common. 

RIT IS TO BE NOTED THAT SIGMA AND RHOL ARE FUNCTIONS OF THE 
RDRO'P' TEMPERATURE , NOT THE VAPOR TEMPlRATURE. IT IS THE DROP 
RTLMPERATURE THEN THAT MUST_ BE USED IN CA LLIN G ON THESE TWO __ 
“RENTRTES, WHEREAS THE VAPOR TEMPERATuRe'Ts USED WHEN 'COMPUT ING 
RTHE LATENT HEAT. 

IT 

EXTERNAL FUNCTION (T) 



PROGRAM common gamma, 


MU* 


B* 


C* CP * 


L, RHOL, SIGMA 






R "■ 

ENTRY TO F RHOL . 














R COMPUTE THE VALUE 

FUNCTION RETURN VALUE 


OF 


THE 


LIQUID 


DENSITY 






R ' 

ENTRY TO FL. 














“R COMPUTE THE VALUE" 

FUNCTION RETURN VALUE 


oT 


The 


LATENT 


HEAT 






R 

ENTRY TO FSIGMA. 














R “'COMPUTE THE" VALUE 

FUNCTION RETURN VALUE 


OF 


THE 


SURFACE 


TENSION 







R 

END OF FUNCTION 



APPENDIX C 



INPUT DATA 



1. Summary of input data and units required. 

2. Table of input data used to obtain results plotted in Figures 1 
through 13. 

3. Plots of vapor properties obtained from References 20 through 
24. 
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Summary of Input Data Requirements 



Parameter 


Definition 


Units 


PZERO 


Initial pressure (P q ) 


2 

dyne/ cm 


RHZERO 


Initial density (p Q ) 


gm/ cm 


TZERO 


Initial temperature (T q ) 


°K 


PE 


Initial approximation to saturation pres- 
sure (P ) 
e 


dyne/ cm^ 


GAMMA 


Ratio of specific heats (y) 


N. D. 


MU 


Molecular weight of the vapor (ju) 


gm/ gmol 


CP 


Specific heat (C ) 


dyne -cm/ gm- 


L 


Latent heat (L) 


dyne-cm/gm 


SIGMA 


Surface tension (cr) 


dyne/ cm 


RHOL 

C b} 


Liquid density (p^) 

Saturation curve constants j 


gm/ cm^ 
dyne/cm^ 
and °K 


ALPHA 


Accommodation coefficient (a) 


N. D. 


DELX 


Incremental step in X (Ax) 


cm 


XRANGE 


Length of interval, starting with conden- 
sation point, over which flow conditions 
will be computed 


cm 


XPOINT 


Points which will be printed out, i. e. , 
every fifth, tenth, etc. , point. 


N.D. 


TRANGE 


Temperature range over which program 
will search for condensation point 


°K 


DELTAT 


Initial step size in condensation point 
search over TRANGE 


°K 


PERCNT 


Value of epsilon for determination of. onset of 
condensation 


Jdecimal 

(fraction 


ASTAR 


Throat area (A*) 


2 

cm 


XMIN 


X- coordinate of intake 


cm 


XMAX 


X-coordinate of exit 


cm 


INANG 


Inlet half angle 


degrees 
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Summary of Input Data Requirements (continued) 



OUTANG 

NA 


Diffuser half angle (6 

Avogadro's number (in program) 
6. 027 x 10 23 


degrees 

molecules/ gmol 


R 


Universal gas constant (in program) 
8. 314 x 10? 


dyne-cm/gmol°K 


K 


Boltzmann’s constant (in program) 
1. 379 xl0“ 16 


dyne-cm/°K 
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N 2 U) 



p 

0 


8. 45 x 10 


p o 


.00952 


T 

0 


298 


p 

e 


2000 


y 


1. 4 


V- 


28. 02 


c 

p 

L 


1. 07 x 10 7 
2. 3 x 10 9 


a 


25. 8-. 22T 


P L 


1. 177-. 00476T 


C 


882 


B 


25.7 


a 


1 . 0 


A* 


. 0645 


X 

mm 


-5 


X 

max 


25 


INANG 


45 



9. 75 



Table C-l. Input Data 


N 2 (II) 


Cu 


8. 45 x 10 6 


1 x 10 6 


. 00961 


2. 54 x 10~ 4 


295 


3000 


2000 


166.7 


1.4 


1. 667 


28. 02 


63. 54 


1. 07 x 10 7 


3. 26 x 10 6 


2. 3 x 10 9 


4.94 x 10 10 


24. 25-. 22T 


1300 


1. 177 ■>. 00476T 


8.9-. 0007T 


882 


34, 200 


25.7 


25. 46 


1 . 0 


1.0 


. 0645 


. 645 


-5 


-5 


25 


58 


45 


45 


9.75 


20 



Zn(I) 


Zn(n) 


1 x 10 6 


1 x 10 6 


3. 93 x 10" 4 


6. 04 x 10' 4 


1500 


1300 


166.7 


166. 7 


1. 667 


1. 667 


65. 37 


65. 37 


3. 18 x 10 6 


3. 18 x 10 6 


1. 82 x 10 10 


1. 82 x 10 10 


1000-. 26T 


1000-. 26T 


7. 4-. 0008T 


7. 4-. 0008T 


13,800 


13,800 


25. 3 


25. 3 


1.0 


1.0 


. 645 


. 645 


-5 


-5 


21. 1 


58 


45 


45 


45 


20 



OUTANG 



P (dynes/cm ) 




Figure C-l. Vapor Pressure of Nitrogen, NBS Circular 564 (24) 



gm/cm dynes/cm 
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Figure C-2. Surface Tension and Liquid Density versus 
Temperature for Nitrogen 
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Figure C-3. Vapor Pressure Curves 



dynes/cm P T gm/cm P T gm/cm 
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°K 



Figure C-4. Metal Vapor Properties 



L cal/g mol cal/g mol 
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Figure C-5. Metal Vapor Properties from Thermochemistry for Steelmaking (20). 



